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ABSTRACT 

We present the dust properties and star-formation histories of local submillimetre- 
selected galaxies, classified by optical morphology. Most of the galaxies are late types 
and very few are early types. The early-type galaxies that are detected contain as 
much dust as typical spirals, and form a unique sample that has been blindly selected 
at submillimetre wavelengths. Additionally, we investigate the properties of the most 
passive, dusty spirals. 

We morphologically classify 1087 galaxies detected in the Herschel-ATLAS Sci- 
ence Demonstration Phase data. Comparing to a control sample of optically selected 
galaxies, we find 5.5% of luminous early-type galaxies are detected in H- ATLAS. The 
H- ATLAS early- type galaxies contain a significant mass of cold dust: the mean dust 
mass is 5.5 x IO^Mq, with individual galaxies ranging from 9 x 10^ — 4 x IO^Mq. This 
is comparable to that of spiral galaxies in our sample, and is an order of magnitude 
more dust than that found for the control early-types, which have a median dust mass 
inferred from stacking of (0.8 — 4.0) x IO^Mq for a cold dust temperature of 25-15 K. 
The early-types detected in H- ATLAS tend to have bluer NUV — r colours, higher spe- 
cific star-formation rates and younger stellar populations than early-types which are 
optically selected, and may be transitioning from the blue cloud to the red sequence. 
We also find that H- ATLAS and control early-types inhabit similar low-density envi- 
ronments. We investigate whether the observed dust in H-ATLAS early-types is from 
evolved stars, or has been acquired from external sources through interactions and 
mergers. We conclude that the dust in H-ATLAS and control ETGs cannot be solely 
from stellar sources, and a large contribution from dust formed in the ISM or external 
sources is required. Alternatively, dust destruction may not be as efficient as predicted. 

We also explore the properties of the most passive spiral galaxies in our sample 
with SSFR < 10~^"'^yr~-'^. We find these passive spirals have lower dust-to-stellar mass 
ratios, higher stellar masses and older stellar population ages than normal spirals. The 
passive spirals inhabit low density environments similar to those of the normal spiral 
galaxies in our sample. This shows that the processes which turn spirals passive do not 
occur solely in the intermediate density environments of group and cluster outskirts. 

Key words: galaxies: classification, colours, luminosities, masses - galaxies: elliptical 
and lenticular - galaxies: evolution - infrared: galaxies - submillimetre - ISM: dust, 
extinction 



1 INTRODUCTION 

It has long been known that there is a relationship be- 
tween galaxy optical colour and morphology. Galaxies can be 
split into a red sequence and blue cloud ([Trcssc ct al. 199^; 
Strateva et alll200ll: iBlanton et ai]|2003l : iBaldrv et al.ll20o1 : 



BelleialJl20o| )~Re d galaxies ar; ire ;ally passive elli7: 



type galaxies (ETGs, those that have elliptical or SO mor- 
phology), but with ~ 25% being spirals which are red eithe r 
due to dust or because they are passive (|Driver et al.ll200^ . 
Blue galaxies are actively star-forming and mostly of spiral 
or irregular morphology. The colour bimodality of galaxies is 
linked to their star- formation history (SFH), with the stellar 
population of galaxies transitioning from blue to red as their 
star-formation ceases due to the r emoval or consumption o f 
cold gas (e.g. iFaber et all l2007l : iHughes fc Cortesd l2009ll . 
The injection of gas and dust via mergers may temporarily 
rejuvenate star- formation, and so this evolution of colour 
can be reversed (ICortese fc Hugh es'2009': ^ Kannappan et al] 
I2OO9I : IWei eralll2010l i. Such "rejuvenators" may have had 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
f E-mail:ppxkr@nottingham. ac.uk 



substantially different star-formation histories from those 
which make up the majority of their morphological type, 
and may provide insight into the evolutionary processes that 
shape galaxies today. 

Massive (> 10^°M b) ETGs are traditionally thought t o 
be "red and dead" fe.g. lTemi. Brighenti fc Mathewsll2009bf) . 
having formed most of their stellar mass a t early epochs 
over a relatively short pe riod of time (e.g. ICimatti et al.l 
I2OO4I : iThomas et al.l |2005| ) and then evolved passively to 
their present state. Their optical light is dominated by old 
stellar populations, however, recent ultra-violet (UV) stud- 
ies of large samples of ETGs have shown that many of these 
galaxies exhibi t low to moderate levels of star-formation 
llYi et al.ll2005l; ISchawinski et al.ll2007l : iKavirai et alll2007l . 
l2008l . I2OIOI : lKavirail[2O10h . The UV-optical colours suggest 
that at least ~30% of UV-selected early-type galaxies at 
z < 0.11 h ave evidence of recen t star- formation within the 
last 1 Gyr l|Kavirai et al.l [20071 ). however, it is difficult to 
determine the contribution of UV flux from old stars. 

Mergers are likely to trigger star-formation, since a high 
incidence of ET Gs with distu rbed morphologies (18%) has 
been observed jKaviraillioTol '). and these disturbed ETGs 
also have bluer NUV — r colours than normal ETGs. The 
major merger rate at low and intermediate redshifts is 
thought to be too low to account for the number of galaxies 
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which have disturbed morp hologies ( e.g. iDe Propris et ahl 
120071 . bOldFLotz et al.ll2o"oi ), therefore IXavirai et al.l (|2010l ) 
conclude that minor mergers are the most likely trigger of 
recent star-formation in ETGs. 

Although there is evidence for limited quantities of dust 
in ETGs, these galaxies are generally thought to be gas- 
and dust-poor, which gives an insight into their evolution- 
ary state. UV starlight is preferentially absorbed and re- 
emitted by dust in the far-infrared (FIR)-submillimetre, so 
the presence of dust emiss ion is often view ed as evidence 
for ongoing star- formation (iKennicuttll 19981 ) . although dust 
can also be heated by the radiation field of an old stel- 
lar population. Ev idence for dust in ETGs was first found 
in th e optical (e.g. Hawarden et al.lll98ll : ISadler fc GerhardI 
1 19851 : Ivan Dokkum fc Fran:>Jll995l ). yet it is difficult to es- 
timate the total dust mass purely from optical observa- 
tions. Warm dust ( > 30K) was detecte d in 12% of lo- 
cal ETGs by IRAS l|Bregman et all 1 19981 ). but IRAS was 
less sensitive to the cold dust co mponent which domi- 
nates the dust mass in local ga l axies ijPunne fc Ealesll200ll : 
IVlahakis. Dunne fc Ealesll2005l : ISmith et al, 1 l2011bl ). There 
have been few studies of ETGs at FIR-submillimetre wave- 
lengths to date, since surveys conducted at these wave- 
lengths have been limited in areal coverage. Consequently, 
studies of ETGs have been targeted observations of rela- 
tively small, and often biased, samples. Cold dust has been 
detected in ETGs throu gh observatio ns with ISO, SCUBA, 
Spitz er and SHARC I I (iTemi et al] e.g. l2004l: iLeeuw et aL 



2004 
2007 



^lahakis e t al.l l2005l: iTemi. Brighenti fc MathewL 
Stickcl, K laas fc Lemkd l2007l : iLeeuw et al.l l2008l : 



Savov. Welch fc FichI l2009i r Cold dust has also been ob- 



served by Herschel in 10 nearby ETGs (jSkibba et al.ll201ll ). 
and in the Virgo cluster elliptical galaxy M86, which con- 
tains dust stripped from the nearby spiral NGC 4438 
ijComez et al.ll201(]l : ICortese et al.ll201(]| ). 

Conversely, spiral galaxies are generally rich in dust 
and gas, and make up the majority of the star-forming 
population. Their blue optical colours indicate young stel- 
lar populations, yet for some time optically red spirals 
with no spectroscopic evidence of star-form ation have been 
known to exist in the outskirts of clusters ([van den Berghl 
1 19761 : iPoggianti et"al]|l999l . l2004l : ICoto et al.ll2003l ). These 
spirals can be red due to dust obscuration , or b ecause 
of an ageing stellar population l|Wolf et al.l |2009| ). It is 
generally believed that passive red spirals have had their 
star-formation quenched due to environmental effects, since 
they are found t o mostly reside in intermediate den- 
sity environments ( Skibba et al.l |2009| : iBamford et al.l l2009l : 
iMasters et~aDl2010bl ). The star-formation rate (SFR) was 
found to be lower for red spirals than blue spirals in 
all environments, which indicates that factors other than 
environment can truncate star-for mation in red spirals 
iBamford et al.ll2009l : [Masters et a"l] r2010b). The same au- 
thors also find that a large fraction of red spirals are massive 
(> 1O'°M0). 

We can now get an unprecedented view of dust in 
local galaxi es from the t he //ersc/ie/- ATLAS survey (H- 
ATLAS. ,Eales et al.ll201ol ). The telescope observes at FIR- 
submillimetre wavelengths across the peak of the dust emis- 
sion, making it an unbiased tracer of the dust mass in galax- 
ies. In this paper we examine the properties of galaxies de- 
tected in the H-ATLAS Science Demonstration Phase (SDP) 



field as a function of morphological type, and highlight inter- 
esting populations which do not conform to the usual trend 
of colour and morphology. In particular, we focus our anal- 
ysis on the properties of H-ATLAS ETGs and how these 
galaxies are different to optically selected ETGs, in addi- 
tion to studying a population of dusty, passive spirals. We 
present the detection of the very dustiest ETGs in a large 
area blind submillimetre survey with Herschel, where the 
lack of pre-selection in other bands makes it the first unbi- 
ased survey for cold dust in ETGs. In Section [2] we describe 
the survey, observations and morphological classifications, 
and present the spectral energy distribution (SED) fitting 
method to explore the properties and SFHs of galaxies in 
our sample; the results of which are presented in Section 
[31 The properties of a population of passive spirals are ex- 
amined in Section |4l The AGN fraction of dusty ETGs is 
explored in Section (5] and we investigate whether environ- 
ment may be an influential factor in the properties of our 
galaxies in Section (6] The submillimetre detected ETGs are 
compared to a control sample of optically selected ETGs in 
Section [T] and we discuss the origin of the dust in Section 
(8] We adopt a cosmology with Q.m = 0.27, Q,a = 0.73 and 
Ho = 71kms-iMpc-\ 



2 OBSERVATIONS AND SAMPLE SELECTION 



The H-ATLAS (|Eales et al.ll201(]| ) is a ~570 d eg^ survey un- 
derta ken by the Herschel Space Observatory l|Pilbratt et al.l 
[ioiO) at 100, 160, 250, 350 and 500Atm to provide an un- 
biased view of the submillimetre universe. Observati ons are 
carried out in parall el mode using the P ACS ( Poglitsc h et al.l 
l201Ci ) and SPIRE ([Griffin et al.i |2010| ) instruments simul- 
taneously. In this paper, we use observations in the SDP 
field, with an area of ~14 deg^ centered on Q=09''05"'30.0^ 
S =00°3 0'00.0" (^ J 2000) . Details of the ma p making can be 
found in [Pascald ([201 ll ): [ibar et al.[ ([2010h . A catalogue of 
^5cr detec tions in any of the 250, 350 and 500/im bands was 
produced ([Rigbv et all [201 1[ ) using the MAD-X algorithm 
(Maddox et al. in prep) and contains 6876 sources. The 5a 
noise levels are 132, 126, 32, 36 and 45mJy per beam at 100, 
160, 250, 350 and 500/im, respectively; the beam sizes are 
~ 9, ~ 13, 18, 25 and 35 arcsec in these bands. 

The H-ATLAS SDP field overlaps with that of the 



Gala x y And Mass As s embly (GAMA) survey ([Driver et al 



'2011; Hill et al.[ [201 1[ : [Robotham et al.l [2010l : [Baldrv et a" 
|2010), which will provide ~ 350 000 spectra for galaxies 
at low redshifts over 6 regions, covering ~ 300 square 
degrees. The GAMA data comprise r-ba nd de fined aper- 
ture matched photom etry as descri b ed in [Hill et al. (20 11f) 
from UV GALEX ([Martin et al.l [2005[ : [Morilssev etll! 
[2007, Seibert et al. in prep.), optical ugriz SDSS DR6 
( Adelman-McCarthv et al.[[200O and near- infrared YJHK 
UKIDSS-LAS ([Lawrence et al.l [2007[ ) imaging. Spectro- 
scopic redshifts and spectra from the AAOmega spectro- 
graph are provided for rpotro < 19.8 or (-K'Kron < 17.6 and 
' modclmag < 20.5) or (zKron < 18.2 and 

''modclmag < 20.50 



J'potro is the r-band Petrosian magnitude, whieh is measured 
using a circular aperture of twice the Petrosian radius, which is 
defined using the Hght profile of the galaxy, rniodolmag is the SDSS 
r-band model magnitude, which is determined from the best fit 
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in the G12 field, and rpctro < 19.4 or (ATKron < 17.6 and 
rmodcimag < 20.5) Or (zKroii < 18.2 and rmodcimag < 20.5) in 
G15 and G09 which includes the H -ATLAS SDP field. 

A likelihood-ratio analysis (|Sutherland fc Saunders! 

Il992f) is performed to match 250/im sources to SDSS DR7 
l Abazaiian et al] 20091) sou rces with r < 22 A within a 10" 
radius ( Smith et al.l 20"llah . and accounts for the possibil- 
ity of the true counterpart being below the optical mag- 
nitude limit. The reliability of an association is defined as 
the probability that an optical source is associated with the 
submillimetre source. SDSS sources with reliability^ 0.8 are 
considered to be likely matches to submillimetre sources, 
these are matched to GAMA survey data to provide spec- 
tra when available. There are 2423 reliable optical counter- 
parts to H-ATLAS sources, with either photometric or spec- 
troscopic redshifts. Around two-thirds of the objects with- 
out reliable optical counterparts are unidentified because 
their counterparts lie below the optical mag nitude limit. 
Thes e sources mostly reside at 2 > 0.5 (see !Dunne et al.! 
!201ll ). The remaining unidentified sources are believed to 
have a counterpart in the SDSS catalogue but we are un- 
able to unambiguously identify the correct counterpart in 
all cases due to near neighbours and the non-negligible 
probability of a background gal axy of the same m agni- 
tude being found at this distance. !Smith et al.l ()2011a! ) esti- 
mate the completeness of the H-ATLAS sample as a func- 
tion of redshift by calculating the total number sources 
that we would expect to have a counterpart above the 
SDSS magnitude limit in H-ATLAS; we refer the reader to 



Smith et al 



Smith et al 



2011al) : [Dunne et al.! (!201l! ) for further details. 
( 2011b ') find that at z < 0.35 the r-band selec- 



tion does not bias our sample towards less obscured sources. 
Since the majority of our spirals and ETGs lie at redshifts 
less than this, our sample should be representative of the 
low-redshi ft galaxy population. Matches ar e also made to 
thelRAS (!Moshir. Kopman fc Conrow "'l992) and FIRST ra- 
di o catalogues (!Becker. White fc Helfandil995i ) as described 
in !Smith et al.! (|2011ah . 



2.1 Morphology 

Morphological classification of sources was performed by eye 
using SDSS standard depth gri composite images, and ob- 
jects were assigned one of four categories: early-type, late- 
type, merger and unknown. The classification fractions are 
shown in Table U] ETGs were identified by looking for a 
dominant bulge and a complete lack of spiral arms, and 
late-types were identified by the presence of spiral arms. 
Due to the shallow depth of the SDSS images, we do not 
discriminate between E and SO types, however, it is pos- 
sible tlmt_th*se_20£ulati^^ properties 
fe.g. !Temi. Brighenti fc Mathew3!2009a! ). The merger cate- 
gory contains systems of galaxies that are clearly interact- 
ing. Galaxies were classified as 'unknown' if it was impossi- 
ble to assign a morphology, usually because the galaxy was 
too faint or small. This situation becomes more common 
as spatial resolution and signal-to-noise decrease at higher 



of an expon ential or de Vaucou leurs profile; further details are 
presented in iBaldrv et al. I ( l20ld) . 



redshifts. It is possible that at low redshifts some of the un- 
known classifications are irregulars, which tend to have small 
angular size and are therefore difficult to identify. Addition- 
ally, very few H-ATLAS galaxies are low stellar mass objects, 
which is due to the flux limit in the submillimetre. Therefore 
the dearth of irregulars is likely to be a real effect and not 
an inability to classify them. Given the sample size, visual 
inspection is the preferred method to classify our galaxies 
into broad morphological classes. It has been shown that 
visual inspection is superior in identifying contaminants in 
samples of ETGs (e.g. face on spirals which have a dominant 
bulge but have weak spiral arms) than automated classifl- 
cation metho ds (Kavirai ct al. 2007; ^hawinski et al.ll2007! : 
!Lintott et ai] I2OO& ; Bamford et al.„ 2009i r Since we are in- 
terested in selecting spheroids, inclination is not an issue. 
It is possible that at higher redshifts Sa type galaxies with 
faint spiral arms not visible in the shallow imaging could be 
classified as ETGs. 



2.1.1 H-ATLAS sample 

We morphologically classify 1087 H-ATLAS sources which 
have reliability ^0.8 of being associated with an SDSS 
source, and which have good quality spectroscopic redshifts 
(flagged with z_quality {nQ) > 3) . Additionally, we require 
that sources are at a redshift oi z < 0.5; above this redshift 
only a very small number of galaxies have spectroscopic red- 
shifts and will be difficult to classify. Two sources with stellar 
or QSO IDs were removed from the sampl e, as were t he fiv e 
sources identified as being lensed in N egrello et al.! (j2010l ). 
We calculate the number of false IDs expected in the clas- 
sified sample from the sum of the probabilities of a false ID 
as 5]](1 — i?), where R is the reliability. This indicates that 
21 galaxies (2%) in our sample are likely to be false IDs. 
There are 115 and 199 sources for which we have PACS 100 
and 160/im point-source detections at ^ 5(7, respectively. 
The selectio n effects aris i ng from the PACS detections are 
discussed in !Smith et al.! (|2011lJ ). who found that the SED 
results in Section l3l are not significantly infiuenced by the 
inclusion of upper limits for PACS data in the majority of 
the sample. All sources are detected ^ 5(7 at 250^m (which 
is a requirement for our sample selection), 272 sources have 
a ^ 5a detection at 350/^m, and 138 sources have a J5 3(t 
detection at 500^m. The distribution of 100/im PACS and 
250/im SPIRE detections is shown for each morphological 
type in Figures IT] (a) and (b). 

We visually classify 44 galaxies as early- type (E or SO), 
with 0.01 < z < 0.32. It can be seen from Table[T]that there 
are few ETGs in our sample compared to spirals, so it is 
evident that H-ATLAS preferentially selects spiral galaxies 
over ETGs. This is as expected since ETGs are generally 
passive and have little dust content. The late-type category 
in principle encompasses both spirals and irregular galaxies, 
however no irregular galaxies are found in our sample. This 
may be because these objects are difficult to classify at all 
but the very lowest redshifts, but H-ATLAS also does not 
detect many low optical lu minosity (and the refore lo w mass) 
sourc es in the SDP field (jDunne et al.![2011 : Darius het al.! 
!201ll ). The number of mergers in our classified sample is un- 
derestimated because the reliability^ 0.8 criteria inherently 
assumes a 1:1 c orrespondence between optica.1 and submil- 
limetre sources (|Sutherland fc SaundersI (l992! : [ Smith et al.! 
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Table 1. Morphologies obtained by visual classification of 1087 H-ATLAS sources and 1052 control sample galaxies. The control sample 
galaxies are selected to have the same r-band magnitude and redshift distribution as those detected in H-ATLAS. The estimated detection 
fraction of galaxies in each morphological class is shown in the last row. These are estimated as explained in Section 12.1.21 
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Figure 1. The distribution of ^ 5cr lOOfim (a) and 250/jm point-source fluxes (b) are shown for the H-ATLAS morphologically classified 
sample. Spirals are shown as grey/filled, ETGs as red/hatched and unknown morphologies as an open histogram. It can be seen that 
the flux distributions are similar for each type of morphology, (c): Redshift distribution of our H-ATLAS sample for each morphological 
type. The ETGs and spirals have similar redshift distributions, but galaxies classified as unknown lie at much higher redshift on average, 
(d): The distribution of SDSS r-band model magnitudes for the same morphological classifications. Sources with unknown morphologies 
are fainter on average than those which are classified. 



l2011al ). In the case of mergers there can be two optical 
sources close to the SPIRE position which both have a high 
likelihood of association but the probability (reliability) is 
split between the sources, sometimes reducing the reliabil- 
ity below our threshold of 0.8. The median redshifts of the 
ETGs and spirals in our sample are both ~ 0.13, and the 
redshift distribution is shown in Figure[T](c). As galaxies be- 



come faint and small with increasing redshift, classification 
becomes difficult, and the unknown fraction increases sig- 
nificantly for Tpctro > 18.5 (see Figure [T] (d)). It also seems 
easier to classify spirals than ETGs at fainter r-magnitudes. 
We observe morphological disturbances in 13/44 (30^5%) 



6 K. Rowlands et al. 



ETGs and 22/496 (4 ± 1%) spiralfl These galaxies show 
evidence of dust or tidal features and may be signs of a 
merger remnant, however this is a lower limit on the num- 
ber which may be disturbed since faint features may not be 
visible in standard depth SDSS images. Morphologically dis- 
turbed sources occupy a range of redshifts up to z ^ 0.26. 
We find a higher f raction of mor phological disturbance in 
ETGs compared to iKavirall (|2010l ) who find 18% for a sam- 
ple of optically selected ETGs (with r < 16.5 and z < 0.05). 

To check our ETG classi fications, we compa r e to t hose 
in the Galaxy Zoo sample (|Lintott et al.1 12008| . HoTJ), in 
which galaxies were visually classified by over 100,000 volun- 
teers. Only the brighter members of our H- ATLAS sample 
(r < 17.77 and z < 0.25) overlap with Galaxy Zoo. Galaxies 
were classified as either ellipticalf], spiral, merger or 'don't 
know'. We assign a galaxy one of these classifications if it has 
> 50 percent of the vote fraction. Debiased votes are used 
to account for the tendency for Galaxy Zoo classifiers to as- 
sign small or faint galaxies (usually at higher redshift due 
to a lack of resolution) to the 'ellipti cal' category. The deb i- 
asing procedure is described fully in lBamford et al"] (|2009t ): 
iLinto tt et al. (201ll). There are 22 of our ETGs which have 
a match in Galaxy Zoo, 17/22 are classified as elliptical, 
3/22 are classified as spiral, and 2 are ambiguous. The ETGs 
which are classified as spirals in Galaxy Zoo either have ev- 
idence of disturbed morphology which could have been mis- 
taken for spiral structure, or have evidence of a disk yet no 
spiral arms. The majority our H- ATLAS ETGs which match 
with the Galaxy Zoo sample are classified as 'elliptical', and 
so our morphological classifications agree well with overlap- 
ping studies. 

We also examine the Sersic index (n) of our ETGs and 
spirals in Figure [2] to check that our morphological classi- 
fications are broadly consistent with what is expected from 
automated galaxy classification. This is accomplished by fit- 
ting single component Sersic models to the light profile of the 
galaxy (Kelvin et al. in prep.). Generally, late-type galaxies 
have an exponential profile {n — 1), and ETGs have a de 
Vaucouleurs profile (n = 4). As expected, our visually clas- 
sified spirals have a very strong peak at n = 1.3, whereas the 
ETGs have a variety of Sersic indicies, but have a higher av- 
erage n of 3.1. The wide range of Sersic indicies is because 
we include SOs in our early-type classification, which may 
have a substantial disk component. We note that although 
n broadly agrees with our visual morphologies, Sersic in- 
dex is not the ideal classification method, because a spiral 
with a bright nucleus may appear to have a high value of n 
and would be mis-classified as an ETG (e.g. iBamford etall 
l2009h . 

2.1.2 Control sample 

In order to understand how the H- ATLAS and optically se- 
lected ETG populations differ, we obtain visual morpho- 
logical classifications of a control sample drawn from the 



0.20 r 




n 



Figure 2. Sersic index (n) distribution of the spirals (grey), and 
ETGs (red/hatched) in our H- ATLAS sample. The dashed line 
at n = 2.5 denotes the traditional cut between 'early-type' (n > 
2.5) and 'late-type' (n < 2.5), and gives an indication of the 
contamination that can occur in samples selected on Sersic index. 

GAMA galaxy catalogue which overlaps with the H- ATLAS 
SDP field. Galaxies are required to be undetected in H- 
ATLAS and have good quality spectroscopic redshifts, and 
were chosen to have the same redshift and rpetro -magnitude 
distribution n(r, z) as our H-ATLAS detected, morphologi- 
cally classified sample. This was accomplished by splitting 
the If- ATLAS sample into (r, z) bins, and randomly pick- 
ing approximately the same number of galaxies in each bin 
from the GAMA catalogue, so that the control sample com- 
prises 1052 galaxies. By selecting a control sample of galaxies 
which are matched in redshift to the H-ATLAS sample, we 
avoid selection effects. 

The morphologies of the control sample are summarised 
in Table [T] It can be seen that there are many more ETGs 
compared to spirals in the optically selected sample than 
in the 250^m selected sample. To estimate the fraction of 
galaxies which are detected at the depth of H-ATLAS as a 
function of morphology, we pick a random sample of 1076 
galaxiefl from the GAMA catalogue in the SDP field, dis- 
regarding whether they are detected by Herschel. We make 
sure the selected galaxies follow the same n{r, z) as our H- 
ATLAS detected and control samples. We repeat the ran- 
dom sampling 1000 times to estimate the average number 
of H-ATLAS detected and undetected galaxies. On average, 
225 galaxies are in the H-ATLAS detected sample, and 11 of 
these are ETGs. There are 851 undetected galaxies, and from 
the control sample fractions we expect 22% (188) of these to 
be ETGs. Consequently, there are 199 ETGs in total in the 
random sample, so we estimate 5.5% of ETGs are detected 
in H-ATLAS compared to the total number of ETGs in the 
SDP field, for this n{r-,z). The detected fractions of other 
morphological types are presented in Table [T] We cannot re- 
liably extrapolate the control sample fractions to the entire 



^ The errors are lu confidence intervals on a binomial population 
using a beta d istribution, whi ch is appropriate for small popula- 
tion numbers || Camero n| l201ll) . 

^ The 'elliptical' cla ssification also c ontain s the majority of SO 
galaxies, as shown in lBamford et ah 1 liooi). 



^ The random sample is chosen to be approximately the same 
size as the H-ATLAS detected sample, but 11 H-ATLAS galaxies 
are not in the GAMA survey region and lack rpotro-magnitude 
information. Therefore the size of the random sample is smaller 
than the H-ATLAS sample, but this should not affect any of our 
conclusions. 
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SDP field, since morphology is a function of both r and z, 
and we have not probed the full (r, z) parameter space in 
this work. 



2.1.3 Classification bias 

iBamford et"al] l|2009l) showed that in Galaxy Zoo the frac- 
tion of galaxies classified as 'elliptical' increases with redshift 
compared to spirals. This is because the spatial resolution 
and signal-to-noise decreases with redshift, so features such 
as spiral arms become invisible. Also, in Galaxy Zoo, images 
are presented to the classifier without any indication of an- 
gular scale, so distant, unresolved galaxies could have been 
classified as elliptical. We should therefore check if this bias 
is present in our classifications. We show the classification 
fractions of our H- ATLAS sample in Figure [S] and there is 
no trend that we classify more ETGs with increasing red- 
shift. Indeed, we see that we classify fewer. This may be 
because unlike Galaxy Zoo volunteers, our expert classifier 
recognises the limitations of the resolution of the image, and 
will classify an object as unknown instead of as an ETG. 

2.2 SED fitting 

ISmith et all (|2011bl ') fit the UV-submilhmetre SEDs of 
1404 H-ATLAS galaxies with reliability > 0.8 of being 
associated with an optical counterpart in the SDSS r- 
band catalogue, and which have available multiwavelength 
photometry. Using the physicall y motivated method of 
IdaTCunha. Chariot fc Elbaj l|2008l . hereafter DCE08) allows 
us to recover the physical properties of these galaxies. In this 
method the energy from UV-optical radiation emitted by 
stellar populations is absorbed by dust, and this is matched 
to that re-radiated in the FIR. Spectral libraries of 25000 
optical models with stochastic star-formation histories, and 
50000 infrared models, are produced at the redshift of each 
galaxy in our sample, containing model parameters and syn- 
thetic photometry from the UV to the submillimetre. The 
optical libraries are produce d using the spec tral evolution 
of stellar populations using a lChabrien l|2003h Galactic-disk 
Initial Mass Function (IMF), calculated from the latest ver- 
sion o f the population synthesis code of lBruzual fc CharlotI 
l|200^ . which includes a revised prescription for thermally- 
pulsing asymptotic giant branch (TP- AGE) stars (Bruzual 
and Chariot, in prep). These libraries contain model spectra 
with a wide range of star-formation histories, metallicities 
and dust attenuation s. The two-component dust model of 
ICharlot fc Fail (|2000l ') is used to calculate the attenuation 
of starlight by dust, which accounts for the increased atten- 
uation of stars in birth clouds compared to old stars in the 
ambient interstellar medium (ISM). The model assumes an- 
gle averaged spectral properties and so does not include any 
spatial or dynamical information. 

The infrared libraries contain SEDs with different tem- 
perature dust components, which include polycyclic aro- 
matic hydrocarbons (PAHs), hot dust (stochastically heated 
small grains, 130 — 250 K), warm dust in birth clouds (30 — 60 
K) and cold dust grains (15 — 25 K) in thermal equilibrium 
in the diffuse ISM, from which the dust mass (Md) is calcu- 
lated. A dust emissivity index /3 = 1.5 is assumed for warm 
dust, and /3 = 2.0 for cold dust, as described in DCE08. 



The attenuated stellar emission and dust emission mod- 
els in the two spectral libraries are combined using a simple 
energy balance argument: that the energy absorbed by dust 
in stellar birth clouds and the diffuse ISM are re-emitted 
by dust in the infrared. In practise, this means that each 
model in the optical library is matched to models in the 
infrared library which have the same value of /^(within a 
tolerance of 0.15) and are scaled to the total dust luminos- 
itjQ Ld°*- We derive statistical constraints on the various 
parameters of the model using the Bayesian approach de- 
scribed in DCE08. We compare each observed galaxy SED 
to the library of stochastic models which encompasses all 
plausible parameter combinations. For each galaxy, we build 
the marginalised likelihood distribution of any physical pa- 
rameter by evaluating how well each model in the library 
can account for the observed properties of the galaxy (by 
computing the goodness of fit). This method ensures 
that possible degeneracies between model parameters are 
included in the final probability density function (PDF) of 
each parameter. The effects of individual wave bands on the 
derived parameters are explored in DCE08, and lSmith et al.l 
(2011b), but we emphasise the importance of using the H- 
ATLAS FIR-submillimetre data to sample the peak of the 
dust emission and the Rayleigh- Jeans slope in order to get 
reliable constraints on the dust mass. 

An example best-fit SED and set of PDFs are shown 
in Figure 2] The parameters we compute are , the frac- 
tion of total dust luminosity contributed by the diffuse ISM; 
TV, total effective V-band optical depth seen by stars in 
birth clouds; Mt/M©, stellar mass; Ld°*/L0, dust luminos- 
ity; Tc^^/K, temperature of the cold diffuse ISM dust com- 
ponent; ty^'^'' , the V-band optical depth in the ambient ISM; 
Md/M0, dust mass; ips/yr^^, specific star-formation rate 
(SSFR); ?/>/M0yr"\ SFR; Ilb, time of last burst; agd, r- 
band light-weighted age and Md/M*, dust to stellar mass 
ratio. For more details of the method we refer the reader to 
DCE08. 



3 PROPERTIES OF ETGS COMPARED TO 
SPIRALS 

Here we explore the multiwavelength properties of our sam- 
ple of morphologically classified spirals and ETGs detected 
in H-ATLAS. We present parameters derived from the SED 
fitting method as described in Section 12.21 for 42 of the 44 
ETGs, and as a comparison we also explore the properties of 
450 out of the 496 spiral galaxies in our sample. We present 
the SDSS images, best-fit SEDs and optical spectra of these 
ETGs in Figure IXT] The galaxies which are excluded from 
our analysis do not have available aperture-matched GAMA 
photometry (2 ETGs, 17 spirals); additionally we reject 29 
galaxies from our analysis which have poor quality SED fits 
with > 30. Our sample covers a range of redshifts, but 
since the median redshifts of the ETGs and spirals are ap- 
proximately the same, differences between the samples due 
to evolution in the redshift range are likely to be small. Ad- 
ditionally, we have checked that the following trends are 
present if we look at galaxies at z < 0.13, and z > 0.13. We 

^ Integrated between 3 — lOOO^m. 
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Figure 3. Morphological classification fractions as a function of redshift for tlie H- ATLAS detected sample (top) and t he control sampl e 
(bottom). Errors bars are the la confidence intervals for a binomial population, derived from a beta distribution fsee ICameron|[201 j) . 
Also shown is the ETG-to-s piral fra c tion f or both samples, where error bars are the Icr confidence intervals for a binomial distribution, 
using the approximation of lGehrel3 l ll98d) . The ETG-to-spiral fraction does not increase with redshift, therefore we do not observe a 
bias towards classifying more ETGs as they become smaller and fainter. 



also observe similar results if we separate our successfully 
classified sample into 'early-type' (n > 2.5) and 'late-type' 
(n < 2.5) using Sersic index. 



3.1 SED parameters 

In order to compare physical parameters for ETGs and spi- 
rals in our sample, we compute the average probability den- 
sity function of parameters derived from our SED fitting. 
The average PDFs of ETGs (red) and spirals (black) are 
shown in Figure O and the mean values and errors are sum- 
marised in Table ICll For each parameter, we use the first 
moment of the average PDF to estimate the mean of the 
population. We can estimate the variance on the popula- 
tion mean as the second moment of the average PDF minus 



the mean squared, divided by the number of galaxies in the 
sample. The error on the mean is simply the square root 
of the population variance. The significance of the differ- 
ence in the means is shown in brackets in Figure [5] and uses 
the quadrature sum of the errors on the mean for the two 
populations. 

The ETGs have a mean of 0.74 ± 0.02, which is sig- 
nificantly higher than that of spirals which have a mean 
ffj, of 0.59 ± 0.01. This means that most of the FIR lumi- 
nosity in ETGs is from dust in the diffuse ISM, which is 
mostly heated by old stellar populations (stars older than 
10^ years). Some ETGs have lower values of indicating 
that more of the FIR luminosity comes from dust in birth 
clouds, which is heated by young stars and implies ongoing 
star-formation. The ETGs in our sample are more massive 
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Figure 4. Top: Example best-fit rest-frame SED of an ETG, with observed pliotometr y points from the UV to the submillimetre (red 
squares). 5<t upper limits are shown as arrows. Errors on the photometry are described in lSmith et~aLl ll2011bD . The black line is the best 
fit model, the green Une is the attenuated optical model, the blue line is the unattenuated optical model, and the red line is the infrared 
model. Bottom: probability density functions of parameters for this ETG. 



than spirals, with ETGs having a mean stellar mass (M*) of 
(4.9lo;g) X 10^*^ Mq compared to M. of (1.9±0.1) x 10^" Mq 
for spirals. ETGs have approximately the same mean dust 
mass (Md) and dust luminosity (L^"') as spirals, although 
the mean ratio of dust to stellar mass (Md/M,) for ETGs 
is lower than that for spirals in our sample by 0.38 dex, 
meaning that ETGs are dust deficient for their stellar mass 
compared to spirals. Our median Md/M, value for ETGs is 
(1.6±0.1) X 10~^, w hich is larger than the average found by 
ISkibba et all (|201ll ) of 1.7 x 10"'' for 10 ETGs. In order to 
contain enough dust to be detected in H-ATLAS, galaxies 
which have low Md/M* need to be more massive in general, 
which may explain why ETGs have a higher M* on average 
than spirals in our sample. We find a mean dust mass of 



(5.5tJ;g) X 10'^ Mq for the sample of ETGs, which is larger 
than the hig hest dust ma s ses fo und in some previous studies 
of ETGs e.g. iTemi et all lj2004h found 10^ - lO'^ Mq, which 
is similar t o the dust m asses found in ETGs with optical 
dust lanes jKavirai et al. 2011) (although these are likely to 
be underestimated by the use of IRAS data - this issue will 
be addresse d in future work . ) Ou r mean dust mass is con- 
sistent with IVlahakis et ahl (|2005l ) who found dust masses 
greater than 10^ M0 for 6 elliptical galaxies from an opti- 
cally selected sample observed with SCUBA. The dust mass 
inferred for the SCUBA ellipticals may include contamina- 
tion from synchrotron radiation, but for the sample of ETGs 
studied here, we find that synchrotron radiation is negligible 
compared to thermal emission from dust (see Section 15. 3p . 
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Figure 5. Average PDFs of the SED parameters of 42 ETGs (red line) compared to 450 spirals (black line). The parameters are (from 
left to right): the fraction of total dust luminosity contributed by the diffuse ISM; M,, stellar mass; M;j, dust mass; M;j/Mt, dust to 
stellar mass ratio; L^°', dust luminosity; Tq^^, temperature of the cold ISM dust component; fy, total effective V-band optical depth 
seen by stars in birth clouds; f^^*^, the V-band optical depth in the ambient ISM; SSFR and SFR averaged over the last 10* years; tLBi 
time of last burst; and agCr, the r-band light- weighted age of the stellar population. The uncertainty on each distribution for ETGs and 
spirals is given by the error on the mean and is shown at the top of each histogram with corresponding colours, and the significance of 
the difference in the means in brackets. The errors for logarithmic parameters are in dex. 



We find no significant difference in tiie Tq of tiie spirals 
and ETGs, and find a wide range of values for tlie dust tem- 
perature. The total effective V-band optical depth seen by 
stars (fv, fy^^') is approximately the same for ETGs and 
spirals. This shows that ETGs have approximately the same 
attenuation as spiral galaxies (though with rather large un- 
certainties). 



3.2 Star-formation histories 

We investigate the SFH of our galaxies by examining the 
SFR, (ip) and SSFR, {tps, defined as tp/M,) averaged over 
the last 10* years. These parameters are derived from the 
SED fitting as described in Section [Ij] The model SFHs 
are described by a continuous exponentially decreasing star- 
formation rate, with sup erimposed randomly dis tributed 
bursts of star-formation (|Kauffmann et al] l2003al ) lasting 
between 3 x 10^ and 3 x 10* years. These bursts occur with 
equal probability throughout the lifetime of the galaxy. The 
probability is set such that 50 percent of the galaxies in the 



library have undergone a burst of star-formation in the last 
2 Gyr. The amplitude of the burst (ratio of mass formed in 
the burst to mass formed in continuous star-formation over 
the lifetime of the galaxy) is distributed logarithmically be- 
tween 0.03 and 4.0. For further details of the models, and 
the effects of mod el assumptions on derived parameters we 
refer the reader to lKauffmann et all |20033) and DCE08. 

The mean SFR for ETGs is 0.7 ± O.2M0yr"\ with a 
range of 0.04—12.4 M0yr~^. It is interesting to note that dis- 
tribution of ETG SFRs in Figure [5] shows signs of bimodal- 
ity. Our range of S FRs is comparable to tha t found for opti- 
cally blue ETGs bv lSchawinski et al.l (|2009l ). who find SFRs 
of 0.5 — 50 M0yr~^ using a range of indicators (Ha luminos- 
ity, u-band light, infrared luminosity from IRAS). H- ATLAS 
ETG SFRs are also larger than those found in recent studies 
of ETGs in the SAURON sample, which is a representative 
sample of lo cal ETGs, which are l ocated in both clu sters 
and the field l|de Zeeuw et al.ll2002l ). lTemi et al.l l|2009ah find 
that the SFR for SAURON SOs as estimated from 24Mm lu- 
minosity is 0.02- 0.2 M0yr"\ and lShapiro et al.l l|2010l ) cal- 



H-ATLAS/GAMA: 



Dusty early-type galaxies and passive spirals 11 



culated the SFR in the SAURON sample from non-stellar 
8/im emission, and this was found to be < O.4M0yr~^. 
These findings of low level star-formation in the SAURON 
galaxies can possibly be explained by the optical selection, 
which is not biased towards highly star-forming galaxies. 
This is in contrast to the H-ATLAS sample which selects 
the dustiest ETGs, and therefore the highest SFRs. Addi- 
tionally, the SAURON measurements only give the obscured 
SFR, and may not be representative of the total SFR of the 
galaxy. 

SSFR is defined as the star-formation rate per unit 
stellar mass and measures the star-formation efhciency of 
a galaxy. Figure [S] shows that the mean SSFR averaged over 
the last 10* years for ETGs {lAt^^l x 10""yr"^) is lower 
than that of spirals (1.0 ± 0.7 x lO'^^yr"!). This trend is 
insensitive to changes in the timescale over which the SSFR 
is averaged. There is, however, a wide range of SSFR and 
17 percent of ETGs have a SSFR greater than the mean of 
the spiral sample. 

In Figure [6] (a), we show a plot of dust mass versus 
SFR for spirals and ETGs in our sample. It can be seen 
that galaxies with the highest dust mass also have a high 
SFR. This trend is expected since both dust mass and SFR 
will depend on the total stellar mass of a galaxy. We can 
remove this trend by dividing by stellar mass and so we 
plot Md/M« vs . SSFR in Figure [6] (b). As was found in 
Ida Cunhaeralll|2010l ). there is a strong correlation between 
these two parameters. It can be seen that typically the ETGs 
have lower SSFR and Md/M* than spirals. There are some 
spirals with very low SSFR and Md/M*, which are discussed 
in Sectional 

We can use the results of our SED fitting to see if 
star-formation is dominated by a recent burst or contin- 
uous star- formation using the model parameter Ilb- Al- 
though there is a large uncertainty on this parameter, our 
results are still useful for a statistical comparison of two 
populations. As shown in Figure [S] ~ 76 percent of our 
ETGs have not had a burst of star-formation in the last 
10^ years, and have therefore not formed a substantial frac- 
tion of their mass in recent bursts. It seems that most of 
our sample have residual star-formation le ft over from the 
last major burst. iKauffmann et all l|2003bl ) find that galax- 
ies with M, > 10^° Mq typically have not had recent bursts 
of star formation, which may explain why our generally high 
mass ETG sample shows few recent bursts. The time since 
the last burst can also be characterised by the age of the 
young stellar population, parametrised in our models by the 
r-band light-weighted age (age^). It is found that the mean 
stellar population age of the ETGs is 2.8 ±0.3 Gyr, which is 
older than that found for the spirals of 1.6 ± 0.1 Gyr. This 
is consistent with the general picture that ETGs are older 
than spirals. We note that 3/10 ETGs with bursts of star- 
formation in the last 1 Gyr show disturbed morphologies, 
so galaxy interactions may be the cause of the burst. It is 
possible that more ETGs in this sample are disturbed at a 
level which is not detected in the shallow imaging that we 
have available. Without deeper imaging, conclusions can- 
not be drawn about whether there is a correlation between 
morphological disturbance and recent star-formation in this 
sample. 
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Figure 6. (a): Md vs. SFR of the ETGs (large red stars) com- 
pared to the spiral galaxies (small grey crosses), (b): Md/M, vs. 
SSFR of the ETGs compared to the spiral galaxies. The error 
bars indicate the median Ic uncertainty on a data point. 



3.3 Comparison of broadband photometric and 
spectroscopic star-formation parameters 

SFH parameters are traditionally measured using spectro- 
scopic information, whereas we used broadband data, so 
there may be a lar ge uncertainty on some parameters. 
IWalcher et al.l l|2008l ) explored de generacies in the SFH pa- 
rameters from broadband photometry using similar stellar 
population models to those in thi s work, and f ound that M« , 
age^ and are well determined. IWild et al.l (|2009l ) classify 
galaxies into star-forming galaxies and quiescent galaxies 
using broadband and spectroscopic data, and found a good 
agreement between these two classification methods. They 
also find the time of last burst derived from broadband SED 
fitting agrees with that derived from spectroscopy. 

To investigate whether fitting SEDs to broadband pho- 
tometry can accurately describe the SFHs of our galaxies, 
we stack spectra of ETGs and spirals together in bins of 
SSFR and r-band weighted age to look for trends in spec- 
tral features. The spectra are shifted to rest wavelength and 
resampled onto a common wavelength array. The spectra are 
normalised to the median of the spectrum, and then com- 
bined using the median of the spectra in each bin. Spectra 
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which show signs of AGN (see Section [5]), or have anoma- 
lous effects such a s bad sky subtraction or fibre fringing 
(jColless et al-lfioOll ) have been removed. It can be seen from 
Figure [7] that as expected, the galax;ies with the highest 
SSFR show signatures of star-formation such as strong Ha 
and [OIII] and [Oil] emission lines. Going from high to low 
SSFR, the strength of the emission lines decrease; the same 
trends are found for age, with older stellar populations show- 
ing minimal signs of star-formation. 



3.4 UV-Optical Colours 

Galaxy colour is often used as a proxy for the age of a stel- 
lar population, with red galaxies assumed to be old due to 
a lack of UV emission from young stars. This simple inter- 
pretation can become complicated, with young star-forming 
galaxies appearing red due to dust obscuration, and old 
galaxies appearing blue due t o con tamination of the UV light 
by horizontal branc h startFI (lO'C onncll 19991: lYi et allbOOSl : 
iKavirai et al.|[2009l ). iDariush et al.. (,2011i ) separate red and 
blue galaxies in the H-ATLAS sample at NUV - r = 4.5 
by fitting double Gaussians to the colour distribution. They 
found that Herschel preferentially selects blue galaxies, and 
that 90 percent of H-ATLAS sources with red colours are 
not old/passiv4l| but have their light attenuated by dust. 
We examine the NUV — r colours of our morphologically 
selected galaxies using aperture matched GALEX UV and 
GAM A optical photometry in Figure [S] Rest-frame photom- 
etry i s calculated using k-correct.v4.2 (|Blanton fc Roweij 
\2m% . and is corrected for galac tic extinction usin g the red- 
dening data of Schlcgcl, Finkbci ner fc Davis! l| 19981 ). Overall, 
93 percent of ETGs have available NUV photometry. For 
sources which have a < 5ct NUV detection 0, we compute 
lower limits for the colours. The mean error in the NUV — r 
colour is 0.08 magnitudes. 

Using the colour cut of lDariush etal\ (|201ll ) at NUV - 
r — 4.5 in Figure [8] (a), we find the ETGs have a range of 
colours, with 24 'blue' and 15 'red' ETGs. Many exist in the 
transition region between the red sequence and blue cloud. 
The SSFR of each galaxy is represented by the colour of each 
point, and a correlation with NUV — r colour is observed. 
As expected, blue galaxies tend to have a higher SSFR, and 
red galaxies a lower SSFR, although with some exceptions. 
In Figure [H] (b) there is a wide range in the colours of both 
morphological types, although the median NUV — r colour 
for the spirals is bluer than that of the ETGs. This trend is 
expected since spirals have the bulk of their stellar popula- 
tion dominated by young stars. 

The red ETGs generally have low SSFR, but still con- 
tain an appreciable amount of dust. These sources have high 
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Figure 8. (a): UV-optical colour magnitude diagram, colour 
coded according to SSFR. Circles are ETGs, crosses are spirals, 
filled circles indicate that the ETG is morphologically disturbed. 
The dash ed line shows the sep aration between 'blue' and'red' as 
defined in toariush et al.l l l201lh . Lower limits are shown for galax- 
ies which do not have a 5o" NUV detection, (b): The distribu- 
tion of (NUV — r) colours for the ETGs (red/hatched) and spirals 
(grey) which have a ^ Scr NUV detection. 



° We note that UV contamination from old stars is unlikely to 
be a con cern, since our sample does not co ntain giant elliptical 
galaxies jYi, Demargue fc Oemler Jr.llT997l) , and UV flux from 
old st ars is likely to be swamped by that produced by young 
st ars llKavirai et al.ll201Cl) . 

' iDariush et al. 1 1I2OIII) define 'passive' systems as galaxies which 
have red {NUV — r > 4.5) colours, after correcting for dust ob- 
scuration. 

* corresponding to NUV > 23.0 in the AB magnitude system 
after galactic extinction correction. 



ffi values which indicate the dust in these sources is predom- 
inantly heated by an older stellar population, which gives 
rise to the red colour of these galaxies. We are observing 
these objects at a time when their star-formation has mostly 
ceased, either because they have used up all their gas, or 
because star-formation has been quenched by some process. 
Their dust has not yet been destroyed by sputtering and 
shocks from type la SNe, and this is discussed in Section [S] 
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Figure 7. Stacked spectra in bins of SSFR (left) and stellar population age (right). The bins are arranged from high SSFR/young stellar 
population (top) to low SSFR/old stellar population (bottom). The grey area shows the normalised, individual spectra, and the thick 
black line is the median of these spectra, smoothed by a boxcar of SA. Prominent emission lines are shown by the blue dashed lines. The 
number of spectra in each stack is indicated in each panel. 



The ETGs which show signs of morphological distur- 
bance (denoted by filled circles in Figure [8] (a)) span a 
range of colours. A Kolmogorov-Smirnov (K-S) test gives 
a probability of 0.14 of the colours of disturbed and non- 
disturbed ETGs being drawn from the same distribution, 
howe ver, this i s no t significantly different (l.lcr). In con- 
trast, iKavirail jioiQ ) find that peculiar ETGs have signifi- 
cantly bluer NUV — r colours than relaxed ETGs. There is 
also a small population of 15 spirals with NUV — r > 4.5, 
and these are discussed in the following section. 



4 PASSIVE AND RED SPIRALS 

There has been much discussion in the literature 
about whether the red colour of some spirals is 
due to dust extinction or an old stellar popula- 
tion (IWolf. Gr av fc Mciscnhcimor' '2005|; IWolf et all |2009| : 
[Masters et al. I [2010b). Wolf ct_al. (,200i) find optically red 
spirals have a lower SFR than blue spirals, but also contain 
large amounts of dust which obscures star-formation. This 
may be due to the inclusion of edge-on spirals in their sam- 
ple, which would inherently have a higher dust extinction 
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because the central dust lane is oriented along our line-of- 
sight. 

Of the 15 red {NUV — r > 4.5) spirals in our sample, 
only two have moderate levels of star-formation with SSFRJS 
lO^^^yr^^. The majority of the red spirals have SSFR much 
lower than this. By selecting spirals with SSFR< 10~^^yr~^ 
we explore the properties of the 19 (~ 5%) most passive 
galaxies in our spiral sample. We not e that this is differen t 
from the 'passive' definition used by iDariush et all (|201ll ). 
which was based on dust-corrected UV-optical colour. The 
error on the SSFR for some passive spirals is large (up to 2.1 
dex), meaning that some passive spirals could plausibly have 
SSFR> 10""yr-\ however, the mean of the average SSFR 
PDF is (2.6to:g) X 10"^ Vr"\ As a population, we can regard 
the average SSFR of passive spirals as being significantly 
(9.3a") different from those of normal spirals (which have 
a mean of (1.2 ± 0.1) x 10"^°yr-^). SDSS images, best-fit 
SEDs and optical spectra of the passive spirals are presented 
in Figure [Bl] 

These spirals have NUV — r colours ranging from 4.3 
to 5.5, although there are 2/19 spirals for which NUV mag- 
nitudes are not measured due to the source being in close 
proximity to a bright star. We find 13/17 passive spirals are 
'red', and 3/17 are 'blue', with one passive spiral having 
ambiguous colour due to an upper limit on the NUV mag- 
nitude. The majority of the passive spirals are not found at 
the extremes of the colour distribution, and lie in the green 
valley. 

4.1 Properties of passive spirals 

A comparison of average PDFs derived from the SED fitting 
for 19 passive and 431 normal (SSFR> lO^^^yr^^) spirals 
is shown in Figure [UJ and the mean values of the parame- 
ters with errors are summarised in Table ICll The passive 
spirals have a high mean of 0.87 ± 0.02, indicating that 
the majority of the dust luminosity is produced in the dif- 
fuse ISM, and powered mostly by old stellar populations. 
The distribution of V-band optical depths in the passive 
and normal spirals is similar, which argues against the pas- 
sive spirals being red due to higher dust obscuration. The 
diff erences found in o pacity between our passive spirals and 
the I Wolf et al.l (|2009l ) red spirals (which have twice the dust 
extinction of blue spirals) may be because we only examine 
passive spirals, and they select their sample of red spirals on 
the basis of optical colour alone. As we have shown in Sec- 
tion [3]4l red colour does not necessarily mean that galaxies 
are passive. 

The mean M« of the passive spirals is (4.2lg g) x 
10^" M©, in comparison to that of the normal spiral pop- 
ulation which has a mean of (1.9 ± 0.1) x 10^'' M©. We 
find 95 percent of the passive spirals are massive with 
M* > 10^" Mq . This could be a selection bias in that we 
can only detect dust in the most massive passive spirals as 
their Md/M* ratios are muc h lower than t h e norm al spi- 
ral population. Alternatively, [Masters et al.l l|2010bl ) found 
that in their sample almost all red spirals were massive 
(M* > 10^° Mq). Figure [9] shows that the passive spirals in 
our sample have much older stellar population s than the nor- 
mal spi ral population. This is consistent with [Masters et aP 
(|2010bl ). who found red, face-on spirals have older stellar 
populations than blue spirals, and are not post-starburst 




-5 -4 -3 -2 -1 8.0 8.5 9.0 9.5 10.0 

logio(Md/M.) Iog,o(oger / yrs) 



Figure 9. Average PDFs of the SED parameters of 19 passive 
spirals (SSFR< 10~^^yr~^) (blue line) compared to 431 normal 
spirals with SSFR^ 10~^^yr~'^ (black line). The parameters are 
(from left to right): /p, the fraction of total dust luminosity con- 
tributed by the diffuse ISM; M*, stellar mass; M^/Mt, dust to 
stellar mass ratio; and ager, the r-band light- weighted age of the 
stellar population. The uncertainty on each distribution for ETGs 
and spirals is given by the error on the mean and is shown at the 
top of each histogram with corresponding colours, and the signif- 
icance of the difference in the means in brackets. The errors for 
logarithmic parameters are in dex. 

objects. This suggests our spirals have not stopped forming 
stars recently, and may have low SSFR because they have 
used up most of their gas. This implies that, under some 
circumstances, spirals can retain their spiral appearance for 
a fe w Gyr following the cessation of their star-formation 
(e.g. iBekki. Couch fc Shioval |2002| ). This interpretation is 
supported by the time of last burst, for which we find a 
mean of I.SIq 4 Gyr for our passive spirals. 

4.2 Inclination effects 

Our sample of passive spirals is separated from the other 
morphologically-classified spirals on the basis of our SED 
fitting results, which uses an 'angle- averaged' approach. 
Results may be biased for sources with high inclinations 
(|da Cunha et al.ll2010l ). so we calculate the inclination of 
our passive spirals to check that there is not a high frac- 
tion of edge-on galaxies in our sample. The minor to major 
observed axis ratio b/a of the SDSS g-band isophote at 25 
mag arcsec"^ can be used to determine inclination. A ratio 
of b/a of ~1 indicates that a galaxy is face-on, b/a decreases 
as the galaxy inclination becomes edge-on. The observed ax- 
ial ratio b/a can be converted into an in clination using the 
relation (as used in [Masters et al.ll2010al ) 

2. (b/g)^ - 
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Table 2. Inclinations (i) in degrees of the 19 passive spirals in 
our sample. The Herschel SDP ID is given in column 1, b/a is the 
minor to major observed axis ratio of the SDSS g-band isophote 
at 25 mag arcsec"-^, foev is an SDSS parameter which is the 
fraction of the galaxy fit by a de Vaucouleurs profile, q is the 
intrinsic axial ratio that would be measured for i = 90°. 
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where q is the intrinsic axial ratio that would be measured 
for an edge-on galaxy (i = 90°). An estimate of q can 
be obtained from the observed distribution of axial ratios 
for SDSS galaxies with different values of the parameter 
foev This SDSS parameter describes the fraction of the 
galaxy light which is fit by a de Vaucouleurs profile (the 
other fraction of the luminosity is fit by an exponential 
profile), and gives information abou t the bulge-to-disk ra - 
tio. We adopt the relation found in [Masters et al.l (|2010al ) 
q = 0.12 -f 0.10 X foev, and use the (?-band defined foev 
The inclinations are listed in Tabled Assuming that galax- 
ies appear approximately edge-on for i > 75°, then a random 
sample of inclinations would lead to 17 percent of galaxies 
appearing edge-on. We find that 5/19 of our passive spi- 
rals have an edge-on inclination, and so within 2a binomial 
errors our sam ple is consistent with a random distribution 
of inclinations. Ida Cunha et al.l (l2010i ) show that the SSFR 
derived from SED fitting may be biased low for high incli- 
nations {b/a < 0.4, corresponding to i > 67°). However, the 
SSFR of these passive spirals are sufficiently low that af- 
ter accounting for this small bias we would still regard the 
majority of these galaxies as being passive. 

We conclude that most of the 'passive' spirals are red 
because they harbour old stellar populations, not because 
of increased amounts of dust w hich obscures star-form ation. 
This agrees with the findings of lMasters et al.l (|2010bl ). who 
find that red spirals have similar dust content (measured 
from Balmer decrements) to blue spirals at the same stellar 
mass. 



5 STAR-FORMATION AND AGN FRACTIONS 
5.1 Emission line diagnostics 

We use optical emission line ratios plotted on a BPT di- 
agram (|Baldwin. PhiUips fc TerlevichllT98ll ) to characterise 
the AGN activity in our H-ATLAS ETGs and spirals. Line 
ratios and equivalent widths (EWs) are derived from the 
SDSS MPA-JH U catalogutFI foemonti et al.ll2004l ) and the 
GAMA survey (|Driver et al.ll201ll '). We regard a line detec- 
tion as > 3ct above the continuum, but lines affected by sky 
emission or fibre fringing are not used. For line fiuxes derived 
from the SDSS sample, corrections are made for stellar con- 
tinuum absorption by subtracting a stellar population model 
from the spectrum, and measu ring emission lines from the 
residual l|Tremonti et al.l |2004| '1 . Where line fluxes are de- 
rived from GAMA measurements, a correction of 1.3 A for 
stellar absorpti on is applied to the EW of the Ha and H/3 
emission lines j Hopkins et al.l 20031: Gunawardhana et al.l 



I2OIII : IWiiesinghe et al.ll201 j T lGunawardhana et all (I2OIII ') 
found for Ha lines with log(Ha EW)< 0.9 there was a dif- 
ference of more than 5% in EW when a range of absorption 
corrections from 0.7— 1.3A was applied. Some of our sources 
are below log(Ha EW)< 0.9, but our results are unchanged 
if this range of absorption corrections are used. In the cases 
where there are multiple measurements of the same galaxy, 
we take the signal-to-noise weighted mean of the line fiuxes. 

We plot the [OIII]/H/3 line ratio as a function of the 
[Nil] /Ha, [SII]/Ha and [01] /Ha line ratios in Figure \W\ 
for spirals and ETGs. We classify objects as AGN or star- 
forming first from the [01] /Ha diagram, since [01] is the 
most sensitive to the presence of an AGN. If the galaxy 
is not present in the [OI] diagram, we use the [SII]/Ha 
diagram, and finally the [Nil] /Ha diagram. On all dia- 
grams, galaxie s that lie abov e the curved line are classified 
as AGN (Kc wlev et al.|[200ll '). and galaxies below the line 
are star-forming. Low signal-to-noise [OIII] and H/? lines 
mean that some sources cannot be located on the BPT di- 
agram. In these cases, a source is classified as an AGN if 
[Nil] /Ha > 0.2. In many of our ETGs we do not detect aU 
four required emission lines, so we derive upper limits and 
locate the galaxy on the BPT diagram where at least two 
lines are presenile. We present the classification fractions of 
ETGs in Table [3l more than half of ETGs are star- forming, 
but 45 percent of ETGs cannot be classified due to their 
weak emission lines. For comparison, spiral galaxies are plot- 
ted on Figure [To] in grey, and lie mostly in the star- forming 
region. 

In a sample of optically selected ETGs, 
'Schawins ki et al.l (|2007l ') found 61% are star-forming, 
and 39% are AGN dominated, which is similar to the 
fractions in our H-ATLAS sample. Since our AGN fraction 
is consistent with that from an optically selected sample, 
this would suggest there is no link between the presence of 
AGN and dust emission, although it is interesting to note 



^ http://www.mpa-garching.mpg.de/SDSS/DR7/ 

In the case where there is H/3 absorption, the H/3 flux is not 
measured in the GAMA spectra, so we derive 3a upper limits. 
Assuming a flat continuum, we estimate the area under a Gaus- 
sian line in pixels (Npix) with FWHM equal to the instrumental 
resolution of 3.5A, and estimate the error on this line given the 
mean noise in the spectrum (ct) as \/Npix X a. 
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Figure 10. BPT diagram showing H-ATLAS spirals (grey) and ETGs (black) with all four lines in each diagram detected at > 3a, with 
error bars shown. Those below the curved line are classified as star-forming, those above the curved line are classified as AGN. Upper 
limits are shown for the ETGs where at least two emission lines are detected. ETGs which have radio emission (Section [STSj are marked 
with a circle, and disturbed morphologies with a diagonal cross. 



Table 3. Emission line classifications of H-ATLAS ETGs which 
can be unambiguously classified on the BPT diagram, with Pois- 
son errors. These fractions do not include galaxies which we can- 
not classify into either category, which comprises 45 percent of the 
sample. Ambiguous classifications result from one or more weak 
emission lines not detected at > 3cr, or measurements affected by 
skylines. The errors are lir confidence intervals on a binomial pop- 
ulation using a beta distri bution, wh ich is appropriate for small 
population numbers (Cameronl l201lh . 



Classification Number ETG Percentage 



ETGs 


23 


100% 


Star-forming 


13 ± 2 




AGN 


10 ± 2 


43t^°% 



that we detect few LINEr£3 

in our sample (although some 
galaxies with upper limits may fall into this category). The 
lack of LINERs in our s ample may be because they are dust 
poor l|Kauffmann et al . 2003b; Kowley et al. 200^, and 
therefore we may potentially be biased against detecting 
LINERS in H-ATLAS, although we need a larger sample of 
galaxies to confirm this. 

We do not account for AGN emission in the SED fitting, 
so we may expect the galaxies with AGN to be poorly fit 
by the models. For the ETGs which host AGN, their SEDs 
look similar in the optical to those w hich are classified as 
star-forming. iKauffmann et ahl (|2003lJ ) find that the optical 
spectra of type-2 AGN have a small fraction of their opti- 
cal light from non-stellar sources, and are very similar to 
spectra of non-AGN host galaxies, except for emission lines. 
Since our physical properties are determined from broad- 
band fitting and not from line strengths, properties from op- 

Low-ionization nuclear emission-line region. 



tical data s hould not be affected by the p resence of a type-2 
AGN. Since iHatziminaoglou et al.l ()201Cl ) find no difference 
between the FIR/submillimetre colours of star- forming and 
AGN galaxies; the FIR is insensitive to presence of AGN 
and therefore will not produce a bias in SED parameters. 

5.2 Hq equivalent widths 

We present the Hq EW distribution of our ETGs in compar- 
ison to spirals in Figure [TT] For the ETGs there is a range 
in EW from 0-109A with a median of 8.7A, which is lower 
than the median for the spirals in our sample (16.4A). The 
median val ue for the spirals is s imilar to that found for field 
galaxies bv lTresse et al.l ()l999( l. It is not unsurprising that 
the EW of ETGs is less than that of spirals, but nonethe- 
less some EWs are substantial and indicate ongoing star- 
formation (consistent with the broad-band SED fitting). The 
ra nge of Hq EWs in ETG s are comparable to those found 
by ISchawinski et al.l (|2009'), who found EWs up to 85A in 
their blue ETG sample. Fukugi ta et al ] (|2004 ) found that 
visually classified ETGs (with r < 15.9 and z < 0.12) have 
a similar Hq EW range as our sample, with 19 out of 420 
E/SOs with Hq EW > lOA (which represents star- forming 
galaxies). In our sample, we find that a much larger fraction 
(31 percent) of our ETGs have Ha EW > lOA, which is 
unsurprising given our submillimetre selection. 

5.3 Radio detections 

Another indicat or of star- formation and AGN activity is 
radio emission. ISmith et al.l l|2011al ) computed the sta- 
tistical probability of a chance alignment between radio 
and H-ATLAS sour ces using the frequentist technique of 
iDownes et al.l (|l986l ) , which used a method to determine the 
most likely radio counterpart by choosing the source with 
the lowest probability P of being a chance alignment. We 
cross-match our ETGs with the FIRST radio catalogue and 
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Figure 11. Ha EW (corrected for stellar absorption) of spiral 
galaxies (grey) and ETGs (red/hatched) which have ^ Scr Ha 
detections. A K-S test shows that the ETGs have a probability 
of 1.7 X 10~* of being drawn from the same distribution as the 
spirals. 



find that 5/42 ETGs have radio counterparts with P < 0.2, 
and so are considered to be likely associations. The radio 
emission may indicate the presence of an AGN and/or star- 
formation, so we compute the ratio of the bolometric IR 
flux to the 1.4 GHz radio flux (pm) u si nR the me thod of 
iHelou. Soifer fc Rowan-RobinsonI (llOSSl VlBelll l|2003i ). qm is 
defined as 



Table 4. qjn values for the 5/42 Herschel ETGs with reliable ra- 
dio counterparts and SED fits. Errors are propagated from the la 
error on TIR and the local noise estimate at the source position 
measured in mjy. Fint is the integrated flux density at 1.4 GHz 
in m Jy. P is the probability of a chance ali gnment of the submi l- 
limetre and the radio source as computed in ISmith et all l l2011al 1. 



Name 


SDP ID 


P 


Fint 


giR 


J091205. 8-1-002656 


SDR 15 


0.014 


4.25 


2.53 ±0.15 


J090352.0-005353 


SDP.45 


0.048 


0.98 


2.59 ±0.17 


J090718.9-005210 


SDP.350 


0.148 


1.18 


2.58 ± 0.18 


J090752.3+012945 


SDP.1027 


0.083 


2.14 


1.68 ± 0.41 


J085947.9-002143 


SDP.6427 


0.077 


8.62 


0.86 ± 0.56 



5.4 Passive spirals 

In most cases the spectra of the passive spirals show little or 
no Ha emission and a strong 4000 A break (see Figure [Blj) . 
indicating low SFR and an old stellar population. Strong 
sodium and magnesium absorption is often observed in the 
spectra, which can indicate the presence of an old stellar 
population, or high metallicity. Only 4/19 passive spirals 
have sufficiently strong emission lines such that they can 
be located on a BPT diagram, and all of these are clas- 
sified as AGN. This may be because AGN are more com- 



mon in massiv e galaxi es (jKauffmann et al.ll2003b ') . although 
[Masters et al. I l|2010bl ) found that red, face-on spirals have 
a higher AGN fraction than blue, face-on spirals. The lack 
of emission lines in the majority of the sample is consis- 
tent with their being selected as passive in terms of star- 
formation, and also indicates a lack of AGN activity. This 
agrees with radio data, as there are no matches for these 
sources in the FIRST radio catalogue. 



giR = logi 



TIR 



3.75 X IQi^Wm- 



logio 



Si. 4 



Wm-2Hz-i 



(2) 

where S'i.4ghz is the rest-frame 1.4 GHz fc-corrected fiux 
density and TIR is the total infrared luminosity (L^°'), 
which is integrated between 3 — 1000/im. 

The qiB. values for the ETGs are presented in Table 
131 Three ETGs have qiR values cons i stent with that found 
for 162 star-forming galaxies in iBeil (J2OO?), with a median 
qjR = 2.64 ± 0.02. We find two ETGs have qm. values which 
are significantly lower than that for star-formation, which 
suggests the presence of a radio-loud AGN in these galax- 
ies. To rule out synchrotron contamination of the 500/im 
flux, we extrapolate the 1.4 GHz radio fiux to 500/im using 
a power law with a spectral slope, a. Assuming a = —0.8 we 
find that the synchrotron emission at this wavelength is neg- 
ligible compared to the dust emission measured at 500/im. 

It is interesting to note the classifications of ETGs us- 
ing emission line ratios are consistent with those from radio 
emission. The three ETGs with radio emission from recent 
star-formation also have some of the bluest optical colours 
and largest Ha EWs, consistent with recent star-formation. 
For the ETGs which are classified as AGN using radio emis- 
sion, one (SDP.6427) is classified as an AGN using emission 
lines. The other (SDP.1027) is likely to be an AGN from its 
line emission, although it has insufficient signal-to-noise to 
confirm this. 



6 ENVIRONMENT OF HERSCHEL 
DETECTED SOURCES 

We examine the environment of ETGs and spirals by com- 
puting the local density around each one (Brough et al. in 
prep). To define the local density, we use a volume limited 
sample of galaxies with Mr < —20 and z < 0.18. The density 
Ejv in Mpc~^ is computed as 



N 



(3) 



where djv is the projected comoving distance to the A'^th 
nearest neighbour within ±1000kms~^, and N = 5. Den- 
sities are computed for all H-ATLAS galaxies which have 
Tpetro ^ 19.4, and have good quality spectroscopic redshifts 
with 0.01 < z < 0.18, which is the limit defined by the 
absolute magnitude limit of the sample. 



6.1 H-ATLAS ETGs 

Using these criteria we are able to measure densities for 30 
ETGs and 354 spirals detected in H-ATLAS, which are com- 
pared in Figure [121 The densities for E TGs and spirals bot h 
range from void to group environments (|Baldrv et al.|[2006l ). 
with most galaxies residing in field environments. There 
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Figure 12. Comparison of densities for H-ATLAS spirals (grey) 
and ETGs (red/hatched). A K-S test (shown at top of histogram) 
shows that we can't statistically rule out the null hypothesis of 
the samples being drawn from the same distribution. 

are few H-ATLAS galaxies in group/cluster environments, 
so our galaxies do not sample the full range of densities 
in the SDP field, which range from (~ 0.01 — 10 0) galax- 
ies M pc~^. From the morphology-density relation (|Dressleil 
1 19801 ). spirals are more numerous in low density environ- 
ments, and ETGs generally reside in high density envi- 
ronments, however, a K-S test reveals there is no signifi- 
cant difference between the densities of spirals and ETGs 
detecte d in H-ATLAS. Th is is consistent with the find- 
ings of iDariush et al.l l|201ll ). who found that the detection 
rate of H-ATLAS galaxies split in to blue and r ed co lours 
does not depend on environment. lYoung et al.l l|20lj ) ob- 
served a volume-limited sample of ETGs and found a sta- 
tistically weak dependence of molecular CO (which is of- 
ten associated with star-formation) on local galaxy den- 
sity, where CO detections were only marginally lower in the 
cluster environment compared to the field. Conversely, blue 
ETGs hav e been found in lower density environments than 
red ETGs (|Schawinski et al. l l2007l : lBamford et al.|[2009l ). al- 
though thes e studies sampled bot h field and cluster envi- 
ronments. [Kann^^a^^iL^D l| 20091 ) found that intermediate 
mass ETGs are common in low density environments, and 
suggest that they may be undergoing disk re-growth. 

It is possible that we are not sampling a large enough 
range of environments to see a significant difference in the 
densities of H-ATLAS sources as a function of morphology. 
The full H-ATLAS data set will encompass the Coma cluster 
and many other rich Abell clusters and will allow a more 
in depth investigation of environmental effects. Since some 
of the H-ATLAS galaxies are in low density regions, it is 
possible that our measure of environment does not always 
reflect the true local density, and instead traces inter-halo 
distances. 

6.2 Passive spirals 

It is thought that passive spirals have had their star- 
formation quenched as a result of galaxy interactions with 



Figure 13. Environment of passive spirals with SSFR< 
10~^^yr~^ detected in H-ATLAS (blue), compared to normal 
spirals with SSFR> lO^'^^yr"^ (grey). A K-S test shows a high 
probability (shown at top of histogram) that both samples are 
drawn from the same distribution. 

the intra-cluster medium. This can remove gas from the 
outer halo, whi ch stops the supply of fuel for star-f ormation 
in the disk (e.g. lBekki eraLll2002l : IWolf et al.ll2009l ). We use 
the local density estimates to test for any environmental dif- 
ferences between our passive and normal spirals, where den- 
sities can be measured for 17/19 passive spirals. Figure [T5] 
shows that passive spirals in our sample mostly inhabit low 
density environments with a median density of 0.28 galaxies 
Mpc~^, which is slightly higher than the median density of 
normal spirals (0.19 galaxies Mpc"'^). A K-S test shows that 
the distributions of densities of passive and normal spirals 
are not significantly diflerent, although this may be due to 
our small sample size . Our median density is different from 
[Masters et al.l (|2010bl) who found the red, f ace-on spiral frac- 
tion peaks at 1 Mpc~^, and lBamford et aL 1 12009) who found 
that the density of red spirals peaks at 6Mpc^^. While 19% 
of their red spirals are found at densities < IMpc"'^, we 
find that 71% of our passive spirals lie at densities lower 
than this. The differences in these fractions may be due to 
selection effects, since the Bamford and Masters samples are 
selected to be 'red', rather than 'passive', and also because 
our H-ATLAS spirals are in low density environments. Our 
sample shows that it is possible to have passive spirals at 
low densities. We can conclude that environment is not the 
only factor influencing whether galaxies are passive, and the 
processes which turn spirals passive occur at both high and 
low densities. 



7 PROPERTIES OF NON-DETECTED ETGS 

We have identified a population of ETGs with substantial 
dust masses, some of which are actively star-forming. In 
order to understand how this population is different from 
optically selected ETGs, we compare to a control sample 
of morphologically classifled galaxies in the SDP field with 
have the same n(r, z) as the H-ATLAS sample, and which 
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are not detected in H- ATLAS. The selection method for this 
sample is described in Section [2. 1.21 

7.1 Dust masses 

Our control sample is comprised of galaxies which are not 
detected in the submillimetre, however, we can investigate 
the average dust mass of optically selected ETGs with stack- 
ing techniques. The stacking is performed on background- 
subtracted, unfiltered SPIRE maps. All detected SPIRE 
sources are subtracted from the map, so we do not contam- 
inate the stack. We stack at the positions of the ETGs in 
the co ntrol sample, using the same method as lBourne et al.l 
l|201ll ). We assume all galaxies are unresolved point sources, 
and for each source we convolve a cut-out of the map with 
a point spread function (PSF) centred on the optical po- 
sition, which is interpolated to the same pixel grid as the 
data map. We share out the flux of blended s ources as de- 
scribed in the appendix of lBouriie et al.ll|201ll ). but find the 
effect to be negligible in this sparsely distributed sample, 
so double-counting of flux does not affect the stacked val- 
ues. This method is effectively similar to stacking in a PSF- 
filtered map. The background level is estimated by stack- 
ing at random positions and this is subtracted from the 
stacked flux. We use the median value in the stack in order 
to avoid bias from outliers. Following the same method as 
iBourne et al.l (|201ll '). we estimate the la error on the median 
from the d i stribu tion of values in the stack as described by 
iGott et al.l l|200ll '). This error estimate automatically takes 
into account both the measurement error, which reduces as 
the square root of the number of objects stacked, and the in- 
trinsic spread of fluxes within the stack. By stacking on the 
positions of 233 ETGs in the control sample, we find me- 
dian fluxes of 2.9±0.5mJy at 250^im (5.8a), 0.8±0.6mJy at 
350/im and — 0.6±0.6mJy at 500/im. The 250/im flux is con- 
si stent with the typica l fluxes of the optically red galaxies 
in lBourne et all (120111 ). 

To obtain the median stacked dust mass, we calculate 
the dust mass of each object in the sample from its measured 
flux in Jy and its redshift, using equationU) Again the error 
is calculated fr om the distr i bution of dust mass values in the 
stack using the lGott et al] l|200ll ) method. 



Md 



S250 DlK 



Ka{v)B{v, Td)(l + 2) 



(4) 



5*250 is the observed 250/im flux, -D_l is the luminosity dis- 
tance at redshift z, B{u, Tk) is the value of the Planck 
function at 250/im for a dust temperature Td, the dust 
mass opacity coe fficient tid(v) is 0.89m^kg~^ (following 
iDunne et al. I2OIII ). K is the fc-correction, which is given 
by 



K ■ 



3+/3 hu^/kTi^, 



ehVa/kTi^a 



(5) 



where is the observed frequency at 250/im, is the emit- 
ted frequency and Tiao is the isothermal temperature of a 
greybody model normalised to recover the stacked flux at 
250/im. We assume a dust emissivity index = 2.0 and 
Tiso = 18.5 K, which ad equately describes t he SEDs of op- 
tically selected galaxies (jBourne et al.ll201ll ). 



Assuming a realistic range of tem peratures of 25-15 kF^ 
dTemi et all |2004| : iLeeuw et all |2004 Smith, M. et al. in 
prep.), we find median dust masses ranging from (0.8 — 4.0) x 
10^ Mq. This is more than an order of magnitude less than 
the dust masses of the H-ATLAS ETGs, indicating that the 
250/im selected ETGs are indeed much dustier than the av- 
erage optically selected ETG. 



7.2 Star formation histories and optical colours 



We use the same technique as described in Section 12.21 to 
fit the multiwavelength SEDs of the control sample galax- 
ies, using 5a upper limits for the FIR-submillimetre fluxes. 
We reject 27 ETGs and spirals which have poor quality SED 
fits with > 30. Although the parameters derived from the 
FIR-submillimetre region of the SED are only constrained by 
the UV-NIR data, we can put similar constraints on SFH 
parameters as for the 250/im selected sample, as most of 
the constraint for SFH parameters comes from the UV-NIR 
photometry. A summary of the parameters derived from the 
mean PDFs is provided in Table ICll Stellar mass is one 
of the main drivers of galaxy properties, so it is important 
to check that the M« distributions are the same for the H- 
ATLAS detected and control ETGs, so that we can compare 
physical properties without a dependence on galaxy mass. 
Figure [141 shows the stacked PDFs of the stellar mass distri- 
butions for H-ATLAS and control ETGs are not significantly 
different, since the control sample is selected to have the 
same r,z distribution. The range of Md/M, for the control 
ETGs is (1.4 - 6.8) x 10"^ for 25-15 K dust, and on average, 
the mean SSFR of the control ETGs is 1.1 dex lower than 
that of H-ATLAS ETGs. A similar trend is found when com- 
paring the mean SFR of the ETGs. For our control ETGs 
the mean r-band light-weighted age of the stellar population 
is 4.6 ± 0.1 Gyr, which is 1.8 Gyr older than the H-ATLAS 
sample of ETGs. 

The NUV — r colours of the control ETGs are com- 
puted as in Section [3.41 and are compared to the H-ATLAS 
ETGs on a colour-magnitude diagram in Figure [15] (a). 
These cover approximately the same range in Mr by de- 
sign. The distribution of colours are shown in Figure [15] (b); 
the control ETGs are on average 1.0 magnitude redder than 
the H-ATLAS detected ETGs. Since the control ETGs are 
not detected in H-ATLAS these galaxies are less obscured 
by dust, with colours dominated by stellar population age 
rather than obscuration. The colour difference between de- 
tected and control ETGs is therefore intrinsic. A handful 
of control ETGs have very blue NUV — r colours, but the 
dust masses of these star-forming galaxies may not be high 
enough to be detected by H-ATLAS. Alternatively, there 
could have been a failure in matching the optical counter- 
part and submillimetre source, which is a possibility for 7 of 
the control ETGs, (of which 3 are 'blue'). These, however. 



Higher dust temperatures have been found in some studies of 
ETGs fe.e. [Savoy et al.ll2009l : ISkibba et al.ll201 j) . but these used 
/3 = 1.5 which results in a higher d ust temperature (~ 3— 4K) be- 
ing calculated l lBendo et aI.|[20o3) . Accounting for this difference 
in /9, these studies yield dust temperatures which are consistent 
with our range of adopted values. 
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Figure 14. Average PDFs of the SED parameters of 42 de- 
tected ETGs (green line) compared to 222 control ETGs (black 
line). The parameters are (from left to right): M«, stellar mass; 
■'Ps/yr'^, SSFR; V/Moyr-^, SFR; and agCr, the r-band light- 
weighted age of the stellar population. The uncertainty on each 
distribution for ETGs and spirals is given by the error on the 
mean and is shown at the top of each histogram with correspond- 
ing colours, and the significance of the difference in the means in 
brackets. The errors for logarithmic parameters are in dex. 



have a very s mall re liability of association as determined in 
ISmith et all l|2011al ). 



7.3 Environments of Herschel non-detected 

sources 

We compare the environments of control sample ETGs and 
spirals in Figure \Wl with densities as calculated in Section 
[6l As expected, on average the median density of control 
ETGs is higher than that of the spirals, and in contrast to 
the H-ATLAS ETGs and spirals, a K-S test shows a low 
probability of the control ETGs and spirals being drawn 
from the same distribution. 

To see how the environments of the detected ETGs are 
different from those in the control sample, we compare the 
densities in Figure [IT] (a), and find they are different at 
only the 1.8a level. There is some indication that H-ATLAS 
ETGs are in lower density environments than optically se- 
lected ETGs, but a larger sample size is needed to confirm 
this. A comparison of the detected and control spirals us- 
ing a K-S test in Figure [T7| (b) shows that we cannot sta- 
tistically rule out the null hypothesis that they are drawn 
from the same distribution. The similarity of the distribu- 
tions suggests that environment does not explain the dif- 
ferences between the H-ATLAS detected and control sam- 
ple ETGs, however, small sample statistics combined with 
a small range of environments currently limits the strength 
of our findings. 
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Figure 15. (a): UV-optical colour magnitude diagram for H- 
ATLAS detected (red stars) and control ETGs (grey crosses). 
Lower limits are shown for galaxies which do not have a ^ 5o" 
NUV detection, (b): Comparison of the NUV — r colours for the 
detected ETGs (red/hatched) and control ETGs (grey /filled) 

which have a ^ 5cr NUV detection. 



8 THE ORIGIN OF DUST OBSERVED IN 
ETGS 

Interstellar dust in early-type galaxies can originate from 
either an internal source ( via mass loss from stars (e.g. 
iGoudfrooii fc de Jond [199a ) or an external source where 
gas and dust are accreted by minor mergers of galax- 
ies or galaxy-gala-xy interactions (e.g . iTemi et al.l l2007l : 



iKavirai et aLll2009l : [siiapiro et al.ll2"oi0r ). Here we use a sim- 
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iMatsuura et al.ll2009l: iDunne et aLlbOHl) and high redshift 




Figure 16. Comparison of densities for control sample ETGs 
(red/hatched) and spirals (grey). As expected, the spirals have 
a lower median density than the ETGs. A K-S test shows a low 
probability of the samples of ETGs and spirals being drawn from 
the same distribution. 



0.20 



0.15 



0.10 



0.05 



0.00 




0.20 r 



0.15 - 



u 0.10 - 



0.05 - 



0.00 



(b) 




-2-10 1 
log,„(S [Mpc-^]) 



-2-10 1 
log,o(E [Mpc-=^]) 



Figure 17. Comparison of environments between H-ATLAS de- 
tected (red/hatched) and control ETGs (a) and H-ATLAS and 
control spirals (b). For the detected and control ETGs a K-S 
test shows the two distributions have a low probability of being 
drawn from the same distribution, although this is only signifi- 
cant at the l.SfT level. We cannot statistically rule out the null 
hypothesis that the control and H-ATLAS spirals are drawn from 
the same distribution. 



pie approach to investigate whether stellar mass loss could 
be responsible for the dust masses of the H-ATLAS ETGs. 

Internal sources of interstellar dust in galaxies are 
poorly understood, particularly the relative contributions 
from supernovae (SNe) and the cool, stellar winds of 
low-intermediate mass stars (LIMS) in their AGB phase. 
In spiral galaxies like the Milky Way (MW), the ma- 
jor d ust source is presumed to be LIMS (e.g. IWhittetl 
120031 '). which inject dust at a rate of 2 x lO"^M0yr"^ 
Dust is destroyed via sputtering in SN driven shocks, with 
theoretical models predicting that dust is d estroyed on 
times c ales of rd e s ^ 30 - 800 Myr in the MW (jjones et all 
1 1994 iTielend Il998l ). Comparing total dust masses 
with destruction timescales for both local (jjonesl I2OOII : 



llMorgan fc Edmunds! I2OO3I: iMichalowski. Watson fc Hiortiil 
I2OIOI : iGall. Andersen fc HiorthI I2OIII) galaxies leads to a 
^dust budget crisis^: the rates of dust injection required to 
maintain the observed dust mass are around an order of 
magnitude higher than the injection rates from LIMS. 

Significant dust production in SNe would alleviate 
this dust budgetary probl em in bot h local and high- 
z galaxies ([Morgan fc Edmulid^ 120031 '). yet there is still 
some controversy about whether SNe can produce the 
~ 1 M0 that would be required to alleviate the bud- 
get problem. Current estimates suggest ^ 10~^ — 1 Mq 



of dust is formed per core-collapse SN llRho et al. 
Dunne et aP 120091: lOtsuka et al.l [20101 : iBarlow et al 



2008 



201c 



Matsuura et al.M201ll . Krause et al. in prep); additionally 
the mass of dust created (if any) in Type la SNe is highly 
uncertain. 

The ISM of ETGs differs from that of spirals due to the 
presence of a hot X-ray emitting component at l(f — 10^ 
K. This presents a harsh environment for dust, and indeed 
the destruction of dust in ETGs is thought to be domi- 
nated by thermal s puttering due to immersion in this hot gas 
l| Jones et al.ll 19941 ). in addition to a contribution from Type 
la supernovae shock collisions. For thermal sputtering in the 
hot ISM, the typical dust destruction timescale is < 50 Myr. 
The destruction timescale is increased in the warm ph ase of 
the ISM (~ 10* K) with Tdcs ~ 1 Gyr l|Barlowl [l973 ). As- 
suming the supernova rate in ETGs is dominated by Type 
la SNe, dust destruction due to Type la shocks depends 
on th e SFR and stellar mass (e.g. IScannapieco fc BildstenI 
l2005h . Estimates for our ETGs (assuming the mean stel- 
lar mass of 4.9 X 10^° M0) are presented in Tabled The 
timescales for destruction from SNe shocks are much longer 
than those from sputtering in hot gas. 

Applying all of these arguments to our sample, (see Ta- 
ble [S] for details) the dust injection rate for an average H- 
ATLAS ETC would need to be 20-400 times higher than 
expected for LIMS and SNe. Taking the full range of de- 
rived values for the SFR for the H-ATLAS ETC sample 
still produces the same discrepancy. This probl em is however 
not co nfined only to our H-ATLAS ETGs, as iDunne et al.l 
l|201lf ) and Gomez et al. (in prep.) found difficulty in pro- 
ducing enough dust in chemical evolution models for all H- 
ATLAS galaxies at 2: < 0.5. 

The median dust mass for the control sample of ETGs 
is approximately 14-69 times lower than the detected sam- 
ple (Md = (0.8 - 4.0) X 10*^ Mq), with a mean SFR of 
0.07 Meyr"^ (Table [5]). If dust is destroyed only by SNe 
shocks then the necessary injection rate is comparable to 
that produced by stars, however if dust sputtering in hot 
gas is important then we still observe more than 10 times 
the amount of dust expected from stars. This suggests either 
sputtering is not as efficient as estimated. X-ray haloes are 
not as abundant or dust reforms quickly in the ISM. It is 
not just the control ETGs which have these issues but also 
the H-ATLAS ETGs. 

In summary, the dust in the ETGs studied here can- 
not have been replen i shed vi a stellar mass loss (see als o 
iGoudfrooii fc de Jond 1 19951 : llrwin. Stil fc Bridged I2OOII ). 
This result is not exclusive to ETGs since this short- 
fall is also seen in local and high redshift galaxies 
l|Morgan fc Edmunds! [20031 ) of all morphologies. Alterna- 
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tively, the dust destruction timescales in ETGs would 
need to be of order 15 Gyrs (two orders of magnitude 
longer than currently estimated). Although this is not 
likely, so me massive ET Gs are seen without a hot X- 
ray ISM l|Pellegrinil Il999l ). and X-ray luminosity is found 
to correlate with the shape of the elliptical light pro- 
file: those ETGs with power law profiles (as opposed 
to cores) a re genera, l ly fou nd to have lower X-ray lu- 
minosities (|Pellegrinil WO^ . Post -merging system s also 



tend to be deficient in hot gas (I Read fc Ponman 199^: 

u7. 



O'S u llivan. Forbes fc Ponman' '200l': 'O'S uUivan fc PonmanI 
2004 : Brassington, Ponman & Read 2007), and a lack of X- 
ray emission was also found in ETGs with stellar popu l ations 
< 4 Gyr old (|Sansom et al.ll2006l ). iBrassington et al.l (|2007l ) 
found that post-merger ETGs within 1 Gyr of coalescence 
are very underluminous in X-rays, and find tentative evi- 
dence that over the period of a few Gyr galaxies regenerate 
their X-ray halo. It may be the case that the dusty ETGs are 
examples of systems with low hot gas content, possibly seen 
within 1 — 2 Gyr after a major star-formation episode which 
created the dust, however, X-ray observations are required 
to confirm this. 

Longer destruction timescales may also be achieved if a 
significant fraction of the star-dust is shielded from the hot 
gas in a cool, dense phase of the ISM. Recent observations 
of ETGs have dete cted significant amounts of cold CO gas 
l|Young et al.|[201ll ). with a detection rate of 22 percent in 
a volume-limited sample of ETGs. It is possible that since 
our ETGs have significant dust mass they also have a large 
cold gas content, however observations of molecular gas are 
required to test this hypothesis. 

Finally, multiple lines of evidence now favour dust 
growth in the cold neutral phase (u h > 30cm~'^) of 
the ISM on r apid timescal es of ~ M yr ('Drainc & Salpeter' 
ll979l : lDwek"fc Scalo 1980; Tielens|[l998 : Drainc 2009). Such 
growth can replenish dust mass lost through shocks and 
thermal sputtering (although seed nuclei must survive sput- 
tering in order for this to occur). Given the shortfall of 
dust from stars, one would then conclude that most (i.e. 
40 - 99% - Table O of the dust in the ISM of Herschel de- 
tected ETGs would need to be grown in the ISM. The dust 
yield of a galaxy would depend only on balancing destruc- 
tion timescales with the formation conditions in the ISM. 

The other scenario we can consider is that mergers with 
dust rich galaxies could account for the dust content of 
the Herschel ETGs . Such a case was demonstrated recently 
when iGomez et al.l (|2010l ) detected ^ 10^ Mq of dust asso- 
ciated with MSG which originated from a recent tidal inter- 
action with the nearby spiral NGC 4438; this is of the same 
order as the dust mass in the stacked, un-detected sam- 
ple. However, to explain the discrepancy between the ob- 
served total dust mass in the _ffersc/ieZ-detected ETGs and 
the dust mass provided via stellar mass loss, we would re- 
quire > 10^ Mq of dust to be accreted. This is akin to the 
average ETG in this sample swallowing a large spiral galaxy 
with an equivalent dust ma ss to the Milky Wa y in its re- 
cent history. As indicated bv lMalleT et all l|2006l '). the major 
merger rate at ^ = 0.3 is 0.018 Gyr"^ (Mgai < 6.4x 10^" Mq) 
to 0.054 Gyr"^ (Mgai > 6.4x 10^° M©) per galaxy. The smaU 
fraction of dusty ETGs (5.5%) relative to similar optically 
selected ETGs could be produced assuming that the dusty 
ETGs are being observed within 0.5 — 1 Gyr after major 



morphological disturbance has subsided l|Lotz et al.ll2008l ). 
on the assumption that sputtering does not destroy the dust 
grains. Observations of gas and stellar kinematics are re- 
quired to test this hypothesis. Lower levels of dust in other 
ETG samples ar e thought to be due to minor mergers (e.g. 
iTemi et al.ll2007D . although in the case of our sample this is 
less plausible since the rate of such mergers needs to be very 
high to create and sustain the dust mass. 

Since only a small fraction of ETGs have dust masses as 
large as those detected in //ersc/iei- ATLAS, it could be that 
our sample represents a short-lived phase in the evolution of 
some ETGs. Either we are seeing them at a time when dust 
is present in a galaxy following a major merger, and the hot 
X-ray component is also suppressed in this phase, allowing 
the dust to survive. Alternatively, we have a sample of ETGs 
with sufficient residual ISM left over from star formation to 
provide a haven for the dust grains to grow and survive the 
hot X-ray gas (which may also be less abundant than average 
in these galaxies). Current dust destruction timescales are 
inconsistent with the amount of dust observed in both LI- 
ATLAS and control ETGs. We can only account for the dust 
observed in some H-ATLAS ETGs if we assume maximal 
supernova dust production (in Type I and II SNe), and no 
destruction, which would produce '-^ 1 x 10^ M© of dust in 
IGyr. 

Further observations to study the kinematics, light pro- 
files, gas content and X-ray properties of this sample will 
be required to fully answer the question of the origin of this 
dust. A more detailed investigation of the evolution of dust 
in these galaxies using a detailed chemical evolution model 
will be presented in Gomez et al (in prep). 



9 CONCLUSIONS 

We present the properties of a 250^m selected sample of 
galaxies according to their morphology. Our sample consists 
of 44 early-type galaxies and 496 spiral galaxies in the 14 
square degree SDP field. Using an energy balance method of 
SED fitting we derive physical parameters, and use these to 
compare the properties of H-ATLAS galaxies as a function 
of morphological type. We also compare to a control sample 
of optical galaxies selected to have the same n(r, z) as the 
H-ATLAS sample. Our main results are as follows. 

• ETGs detected by Herschel are atypical compared to 
optically selected ETGs. We detect significant dust masses 
in H-ATLAS ETGs, with a mean of 5.5 x lO'^M©. Through 
stacking we find that dust masses are an order of magnitude 
lower in optically selected ETGs of a similar stellar mass. 

• Only a small fraction of H-ATLAS ETGs (24 percent) 
have evidence for a recent burst of star formation within 
the past Gyr. Some of these galaxies may have had star- 
formation triggered as a result of an interaction, indicated 
by disturbed morphologies in 31 percent of the sample, al- 
though not all disturbed sources show signs of a recent burst. 
The majority have residual low-level star-formation left over 
from the last burst a few Gyrs ago, and their optical colours 
suggest they exist in the transition region between the blue 
cloud and the red sequence. 

• We find that the control ETGs have lower SSFRs and 
older stellar population ages than H-ATLAS ETGs, which 
is consistent with the red UV-optical colours of the control 
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Milky Way 


LMC^ 


ETGs (H-ATLAS) 


ETGs (control*) 


Average observed dust mass (Mq) 


~ 3 X 10^ 


2 X lO'' 


5.5 X 10^ 


(0.8- 


- 4.0) X 10^ 


Average star formation rate (Mq yr" -'^) 


~ 1.0 


0.25 


0.7 




0.07 




Destruction via SNe T^j^g (Gyr) 


0.5 


2.5 


1.2 




5.0 




Median dust injection rate ^^dl'^dcs (^oy^"^) 


0.06 


8 X 10^"' 


0.045 


(1.6- 


- 8.0) X 10-" 


Range observed dust masses (Mq) 






(3.1 


- 15.0) X 10'^ 






Range star formation rate (Mq yr" ^) 






0.04 


- 12.4 






Range destruction via SNe T^gs (Gyr) 






21 - 


0.07 






Range dust injection rate Mj/tjcs (^Qyr" -"-) 






1.5 X 10~* - 2.2 






Injection rate AGB stars M^^ agb (Mq yr"^)'^ 


2 X 10-3 


5 X 10-5 


(1- 


500) X 10-* 


1.4 X 


10-" 


Injection rate SNe Afin.SNc (MQyr"^)'' 


1 X 10-3 


5 X 10^5 


(0.5 


- 250) X 10-'' 


0.7 X 


10-" 


Time to build up dust mass (Gyrs)*^ 


10 


19 


200 - 


- 2 






Time to build up average dust mass (Gyrs)'' 






400 - 


- 1 


3.8 - 


19.1 


Contribution of dust mass from stellar sources 


5% 


13% 


0.3 - 


-60% 


100 - 


17% 


Extra dust per SNe needed (Mq)^ 


1.0 


0.6 - 0.9 


10 - 


0.1 


0.0 - 


1.01 



Table 5. Summary of the dust sources an d parameters in the H-ATLAS ETG sample with the Milky Way and Large Magellan ic Cloud 
shown for comparison. - llMatsuura et a l. 2009); ^ median dust maa s of control sample for range of temperatures (15-25 K). - ijWhitteti 
I2OO3I . and references therein): - llRho et al.ii200S : lBarlow et aLll2010h for core-collapse SNe; " - calculated from M(j/(Min_AGB + Min^SNe) 
assuming no destruction; ^ - assuming SNe are producing the dust required to match the dust injection rate. 



ETGs. It is possible that the dust content may therefore be 
related to the time of the last major star formation episode 
several Gyrs ago. 

• No significant difference is found in the environments 
of H-ATLAS and control ETGs, although this may be due 
to small sample size. Environment does not seem to influ- 
ence whether an ETG is dusty, at the moderate-low densi- 
ties probed in this study. Additionally, we do not find any 
H-ATLAS ETGs in high density environments. 

• The 5.5% of ETGs detected in H-ATLAS (compared 
to an optical sample of similar stellar mass) contain more 
dust than can be accounted for by production in AGB stars, 
although this problem also extends to the control ETGs (ex- 
cept at the very lowest median dust mass). Most of the dust 
must be formed in the ISM, or an external source of dust 
from major mergers is also a possibility. It is also possible 
that in H-ATLAS and control ETGs the dust destruction 
timescale is longer if they are deficient in X-ray gas. Future 
studies of the kinematics of the gas and stars would be ben- 
eficial in testing this hypothesis, in addition to X-ray and 
CO observations. 

• We examine the properties of passive spirals in our sam- 
ple which have low SSFR< 10~'^'^yr-^, but still contain sig- 
nificant dust mass. They have larger M, and lower Md/M* 
than 'normal' spirals, and are red in colour, which is due to 
an old stellar population, and not due to increased dust red- 
dening. It is possible that these passive spirals have simply 
run out of gas to fuel star-formation, or their star-formation 
has been quenched by some process in the low density envi- 
ronment in which they reside. 



through rejuvenated star-formation, or through the cessa- 
tion of star-formation as the supply of gas ends. 
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Using Herschel we can probe the dust content of dif- 
ferent types of galaxies over a wide range of redshifts. The 
full coverage of the H-ATLAS survey will allow more inves- 
tigation, with far larger numbers of ETGs and passive spi- 
rals. This will improve our understanding of objects which 
are transitioning between the blue and red sequence, either 
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APPENDIX A: EARLY-TYPE GALAXIES 



Table Al. Properties of ETGs derived from SED fitting. The columns are (from left to right): ID, SDP ID, redshift, SDSS RA, SDSS DEC, 250Atm flux in Jy, the fraction of total 
dust luminosity contributed by the diffuse ISM; fy, total effective V-band optical depth seen by stars in birth clouds; Mt/M©, log(stellar mass); L^°*/Lq, log(dust luminosity); Tq^^/K, 
temperature of the cold ISM dust component; flp'^ , the V-band optical depth in the ambient ISM. Mj/Mq, log(dust mass); i/ig/j/r"^, log(SSFR); i/i/MQyr— 1, log(SFR), tLBi log(time 



of last burst); agCr, log(r-band light- weighted age of the stellar population), rest-frame NUV-r colour (Section 13.411 , density (S/galaxies Mpc 
stellar absorption of 1.3A if > 3cr detection). * indicates morphological disturbance. 
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Figure Al. Optical images, multiwavelcngth SEDs and optical spectra of the 42 ETGs in our sample. Images are 30" on a side. The 
rest-frame SEDs of each ETG ar e shown, whe r e red p oints are the observed photometry, with 5cr upper limits shown as arrows. Errors on 
the photometry are described in ISmith et al] l l2011bl) . The black line is the total best fit model, the green line is the attenuated optical 
model, the blue line is the unattenuated optical model, the red line is the infrared model. Spectra are from SDSS and GAMA, and the 
standard deviation in the spectra is also shown. The spectra have been smoothed by a boxcar of 8 pixels. 
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Figure Al - continued 
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Figure Al - continued 
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Figure Al - continued . The wave-like features in the spectrum of SDP.3834 are due to fibre fringing l|CoIIess et al.ll200ll 'l. 
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APPENDIX B: PASSIVE SPIRALS 



Table Bl. Properties of passive spirals derived from SED fitting. The columns are (from left to right): ID, SDP ID, redshift, SDSS RA, SDSS DEC, 250Atm flux in Jy, the 
fraction of total dust luminosity contributed by the diffuse ISM; Ty, total effective V-band optical depth seen by stars in birth clouds; M»/Mq, log(stellar mass); h^°^/hQ, log(dust 
luminosity); T^J^^/K, temperature of the cold ISM dust component; Ty^^' , the V-band optical depth in the ambient ISM. M^/Mq, log(dust mass); ips/yr~^, log(SSFR); ip/MQ-yr—l, 
log(SFR), tLB, log(time of last burst); agei-, log(r-band light-weighted age of the stellar population), rest-frame NUV-r colour (Section l3.4l l. density (S/galaxies Mpc"'^, see Section[6ll 
Ha EW/ A (corrected for stellar absorption of 1.3 A if > 3cr detection). 
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Figure Bl. Optical images, mult i wavelength SEDs and optical spectra of the 19 passive spirals in our sample. Images are 40" on a side. 
The rcst-framc SEDs of each passive spiral arc sh own, where r ed poin ts are the observed photometry, with 5cr upper limits shown as 
arrows. Errors on the photometry are described in lSmith et al.| | |2011b| ). The black line is the total best fit SED model, the green line is 
the attenuated optical model, the blue line is the unattenuated optical model, the red line is the infrared model. Spectra are from SDSS 
and GAMA, and the standard deviation in the spectra is also shown. The spectra have been smoothed by a boxcar of 8 pixels. 
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Figure Bl - continued 
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Figure Bl - continued . The wave-like features in tiie spectrum of SDP.4639 are due to fibre fringing ijColless et al.ll200lh . 



H- ATLAS /GAMA: Dusty early-type galaxies and passive spirals 



APPENDIX C: SUMMARY OF MEAN 
PHYSICAL PROPERTIES 



40 K. Rowlands et al. 



Table CI. Summary of mean physical properties derived from stacking of PDFs for the different populations studied in this paper. The 
parameters are: f/j,, the fraction of total dust luminosity contributed by the diffuse ISM; M^/Mq, log(stellar mass); Mj/Mq, log(dust 
mass); Mj/M*, log(dust to stellar mass ratio); Lj°*/L0, log(dust luminosity); T^^^/K, temperature of the cold ISM dust component; 
fv, total effective V-band optical depth seen by stars in birth clouds; fy^'^ , the V-band optical depth in the ambient ISM, V's/y"^, 
log(SSFR); ip/MQyT~-^ , log(SFR), tLBi log{timc of last burst); agCr, log(r-band light-weighted age of the stellar population). For each 
parameter, we use the first moment of the average PDF to estimate the mean for the population. We can estimate the variance on the 
population mean as the second moment of the average PDF minus the mean squared, divided by the number of galaxies in the sample. 
The error on the mean is simply the square root of the population variance. The errors for logarithmic parameters are in dex. The mean 
parameters from the infrared part of the SED and energy balance parameters are not determined for the control sample, since we only 
have constraints from upper limits on the FIR-submillimetre flux. 



Parameter 


H-ATLAS spiral 


H-ATLAS ETG 


Normal spiral 


Passive spiral 


Control spiral 


Control ETG 




0.59 ±0.01 


0.74 ± 0.02 


0.58 ±0.01 


0.87 ± 0.02 






M. 


10.29 lb 0.02 


10.69 ±0.08 


10.27 ±0.02 


10.62 ± 0.07 


10.15 ±0.03 


10.77 ±0.03 


Md 


7.72 ± 0.02 


7.74 ± 0.08 


7.73 ± 0.02 


7.47 ± 0.10 






Md/M. 


-2.57 ± 0.02 


-2.95 ±0.07 


-2.54 ±0.02 


-3.16 ± 0.09 






T tot 


10.53 ± 0.02 


10.48 ±0.07 


10.55 ±0.02 


10.14 ±0.09 






T-.ISM 

-■c 


19.7 ±0.1 


19.8 ±0.5 


19.7 ±0.1 


19.8 ±0.6 






fv 


2.28 ±0.07 


2.28 ± 0.23 


2.28 ±0.07 


2.34 ± 0.37 


1.66 ± 0.08 


1.61 ±0.10 


ti-ISM 


0.47 ±0.01 


0.43 ± 0.04 


0.48 ± 0.02 


0.41 ± 0.05 


0.24 ± 0.01 


0.20 ±0.01 


V's 


-9.99 ±0.03 


-10.85 ±0.14 


-9.92 ± 0.03 


-11.59 ±0.18 


-10.58 ±0.07 


-11.92 ±0.07 


v> 


0.30 ±0.03 


-0.16 ±0.12 


0.36 ± 0.03 


-0.97 ±0.19 


-0.43 ±0.05 


-1.16 ± 0.07 


tLB 


8.70 ±0.07 


9.04 ±0.18 


8.68 ± 0.08 


9.26 ±0.10 


8.87 ±0.07 


9.39 ± 0.03 


ager 


9.21 ±0.02 


9.45 ± 0.05 


9.19 ±0.02 


9.59 ± 0.05 


9.32 ± 0.02 


9.67 ±0.01 



